March 1995 



CBPF-NF-010/95 



in 



(N 



> 
On 
On 
CO 



5 



THE POMERON: 

YESTERDAY, TODAY and TOMORROW t 



EUGENE LEVIN*) 

LAFEX, Centro Brasileiro de Pesquisas Fisicas ( CNPq) 
Rua Dr. Xavier Sigaud 150, 22290 - 180 Rio de Janeiro, RJ, BRASIL 

and 

Theory Department, Petersburg Nuclear Physics Institute 
188350, Gatchma, St. Petersburg, RUSSIA 
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(i) the brief review of the Reggeon Calculus, which was the first attempt to build the 
effective theory of the strong interaction at high energy. In spite of the fact that this 
approach turns out to be inconsistent and in lectures we show why, the Reggeon approach 
was and is the main source of the terminology and phenomenology for high energy "soft" 
interactions. 

(ii) the detail description of the QCD approach to high energy interaction. We try 
to combine the rigorous approach in perturbative QCD with more simple, intuitive guess 
based on general properties of QCD to clarify our expectations and predictions. 

(iii) the outline of my personal opinion what problems will be important in the future. 

The main reason for the lectures was just the last part to agitate you to think about 
this difficult but because of that interesting problem. The motto of my lectures is " there 
is nothing more exciting than to solve a difficult problem. Pomeron is the one. " 
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Introduction 



In these lectures I want to give the brief review of the main ideas and approaches to 
the structure of the Pomeron. For young generation even the word "Pomeron" sounds 
not familiar. It makes very difficult to explain in the introduction what we are going to 
discuss. In some sense it is the reason why we think that we shall start from the short 
review of the situation with the strong interaction in the past because the terminology 
that we still use came from so called Regge approach which was in the past the common 
language for the phenomenology of high energy strong interaction. 

This is why instead of the usual introduction in which I have to explain the subject 
of the talk I am going to provide you some general outline of the main ideas of the talk. 

I'll start from the review of the Regge approach and I'll try to give you the basics of 
this approach, namely the theoretical background and the main phenomenological ideas 
that had been used in this approach. Just from the beginning I would like to stress that 
inspite of the pure phenomenological input the Reggeon approach was and is the main 
source of the ideas even now for the "soft" high energy interaction. The main reason for 
this is the rich input from the understanding of the general properties of the scattering 
amplitude based on the analyticity and symmetry of strong interaction. The analysis of 
the theoretical problems made at that time ( more than 30 years ago) is still fresh and 
the principle problems which have been pointed out are still with us. 

After the introduction in the past and present problems with the "soft" interaction I 
am going to present what is known about high energy interaction in perturbative QCD 
which is the microscopic theory now. We will go back to old problems and I'll show you 
what kind of understanding we have reached in framework of perturbative QCD (pQCD). 
This understanding consists of two things: the sructure of the Pomeron in perturbative 
QCD ( so called the BFKL Pomeron) and the theory of shadowing correction in the deeply 
inelastic processes. The main idea of this part of my lectures is to show you that we are 
on the right road but only in the beginning. The principle problems are non perturbative 
ones and we have made only the first try to attack them. 

The last, third part will be devoted to outline of my personal opinion what problems 
will be important in the future. You can easily guess that the main reason for these lectures 
was just to agitate you to think about these difficult but really interesting problems. 

As far as concerned the style of these lectures, we tried to follow the Landau ten 
commandments of theoretical physics, at least, the first three of them: 

• Only prediction is a theory, only calculation is a result. Don't believe in a qualitative 
"theory" . 
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• Try to solve a problem exactly, if not try to find a small parameter, if not try again. 

• A model is the theory which we apply to the kinematic region where we cannot 
prove that the theory is wrong. Don't be afraid to build a model, this is the only way to 
be original. 

We would like to make several comments on the quotations. In the first part of the 
lectures we avoid the references except the most important ones. The reason is twofold: 
first is the fact that everybody can find detail description of the Reggeon approach in 
any textbook written in 70's. The second is the idea to give you a full information in the 
body of the lectures about all result that we are going to take with us in the future. In 
the second part we give sufficiently full list of references. 

Concluding this introductory section we want to remind you the motto of these lec- 
tures: 

There is nothing more exciting than to solve a difficult prob- 
lem. Pomeron is the one. 
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YESTERDAY: 

REGGEON APPROACH 

1 The basic ideas. 

1.1 Our strategy. 

To understand our theoretical background in the past let me remind you that before the 
year 1974 the study of the high energy asymptotic was the high priority job, because we 
believed that: 

Analyticity + Asymptotic = Theory of Everything. 

Roughly speaking we needed asymptotic at high energy to specify how many substructions 
we have to make in the dispersion relations to calculate the scattering amplitude using 
them. 

Now situation is quite different, we have good microscopic theory ( QCD) and certainly 
we have a lot of problems in QCD which have to be solved. High energy asymptotic is 
only one of many. I think it is time to ask yourselves why we spend our time and brain 
trying nevertheless to find the high energy asymptotic in QCD. My lectures will be an 
answer to this question, but I would like to start recalling you the main theorems for the 
high energy behaviour of the scattering amplitude which follow directly from the general 
property of analyticity and crossing symmetry and should be fulfilled in any microscopic 
theory including QCD. 

1.2 The great theorems. 

1. Optical Theorem 

Our amplitude is normalized so that 

^ = -|/(M)| 2 (1) 
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Within this normalization the optical theorem says: 

<r tot =AnImf(s,0) (2) 

Therefore the optical theorem gives us the relationship between the behaviour of the 
imaginary part of the scattering amplitude at zero scattering angle and the total cross 
section that can be measured experimentally. 

2. The Froissart boundary. 

We call the Froissart boundary the following limit of the energy growth of the total 
cross section: 

o-tot < C ln 2 s (3) 

where s is the total energy of our elastic reaction: a(p a ) + b(pt>) — > a + b, namely s = 
{Pa + Pb) 2 - Coefficient C has been calculated but we do not need to know its exact value. 
What is really important is the fact that C oc p, where k t is the minimal transverse 
momentum which can be in our reaction. Since the minimal mass is the mass of pion in 
the hadron spectrum the Froissart theorem says that C « ^. The exact calculation 
gives C = 60mbn. 

3. The Pomeranchuk Theorem. 

The Pomeranchuk theorem is the manifestation of the crossing symmetry, which 
can be formulated as a following statement: one analytic function of two variables s and 
t describes the scattering amplitude of two different reactions a + b ^a + baks>0 and 
t < as well as a + b — > a + b at s < 0(u = (p a + Pb) 2 > and t < 0. 

The Pomeranchuk theorem says that the total cross sections of the above two reactions 
should be equal to each other at high energy if the real part of the amplitude is smaller 
than the imaginary its part. 

• 

Problem 1: Prove the Pomeranchuk theorem, using the dispersion relation for the 
scattering amplitude at t — 0. 

1.3 Unitarity constraint in impact parameters (b t ). 

The scattering amplitude in 6 t -space is defined as 

a(s,bt) = ± J d 2 q e-^f(s,t) (4) 
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where t = —q 2 . In this representation 

atot = 2 J d\ Ima{s,b t ) (5) 

a el = Jd\ \a(s,b t )\ 2 (6) 

Using the above notations we can write the s- channel unitarity in the general form: 

2Ima(s,b t ) = \a(s,b t )\ 2 + G in (s,b t ) (7) 

where a(s, b t ) is elastic amplitude while G in is the contribution of the all inelastic processes. 

It should be stressed that the above constraint has a general solution if we assume 
that the scattering amplitude is pure imaginary at high energy, namely 

a(s,b t ) = t(l - e~ n{s ' bt) ) (8) 

where the opacity Q(s, b) is a real function and has a very simple physical meaning. 
Indeed, 

G m {sA) = 1 - e- 2n ^ . 
Therefore e~ 2n is the probability that the incoming particle has no inelastic interaction. 



Problem 2: Prove the Froissart theorem, using the general solution of the s - channel 
unitarity taking into account that at large values of b t the scattering amplitude is small 
and behaves as s N e~^ bt , where // is the mass of the lightest hadron. 

1.4 The first puzzle. 

The first puzzle can be formulated as a question: What happens with resonances with the 
values of their spin bigger than 1 in exchange?" On one hand such resonances have been 
observed experimentally, on the other hand the exchange of the resonances with spin j 
lead to the scattering amplitudes which are proportional to s- 7 where s is the energy of 
two colliding hadrons. Such a behaviour contradicts the Froissart boundary. It means 
that we have to find the theoretical solution of this problem. 

Let me illustrate the problem considering the exchange of vector particle ( j = 1) ( 
see Fig.l). 
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Introducing Sudakov variable we can expand each vector in the following form: 

q = a q p[ + (3 q p' 2 + q t 

where 

p? = P'i = o 

and 

pi = Pi + Pp'2 ; pi = a 2Pi + p'2 ; 

It is easy to find that 

pi — — ; — — 



since 

Pi = m i = 2 PiPi ■ Pi = Pi s 

From equations 

P'l = Ml ■ p'i = M| 
we can find the values of a q and j3 q . Indeed, 

Pi = (p q - qf = ( (1 - a q ) P [ + (Pi - (3 q )p' 2 f = (1 - a q )(Pi - (3 q )s - q 2 t = M 
and 

p'i = (l + P q )(a 2 + a q )s - q 2 t = M\ 

Therefore 

Ml -m\ m\- M\ 

Oi q , p q 

s s 
It is easy to get from the above equations that 

i 2i i n 2i (M% - ml)(Ml - m\) 2 2 
|g 2 | = |a^s-g t 2 | = ^ ^ + g 2 -> |™ q\ 

s 

Now we are prepared to write to expression for the diagram of Fig.l: 
where my is the mass of the vector meson. 



Problem 3: Using the above example show that the exchange of the resonance with 
spin j gives the scattering amplitude equals 

(4sV 

A(s, t = -q 2 t ) = g x g 2 - T - — T 

qf + m 2 R 



Problem 4: Show that the exchange of the resonance with spin j leads to the ampli- 
tude in bt which behaves as s J exp(—mnbt) at large values of b t and s. 



1.5 Reggeons - solutions to the first puzzle. 

The solution to the first puzzle has been found. It turns out that the exchange of all reso- 
nances can be described as an exchange of the new object - Reggeon and its contribution 
to the scattering amplitude is given by the simple function ( see Fig. 2): 

s <*(*) ± (— s )a(t) 

A R (s,t) = g 1 (m 1 ,M 1 ,t)g 2 (m 2 ,M 2 ,t) ' (13) 

sinira{t) 

a(t) is a function of the momentum transfer which we call the Reggeon trajectory. The 
name of the new object as well as the form of the amplitude came from the analysis of the 
properties of the scattering amplitude in t channel using angular momentum representa- 
tion. However, today it is not important the whole history of the approach. What we 
need to understand and take for the future studies are the main properties of the above 
function which plays the crucial role in the theory and phenomenology of the high energy 
interaction. 



2 The main properties of the Reggeon exchange. 
2.1 Analyticity. 

It is obvious that the Reggeon exchange is the analytic function in s, which has the 
imaginary part equals to 

± 9l g 2 s a ® 
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in the s - channel and the imaginary part 



9i 92 s a(t) 

in the u - channel ( for a + b — > a + b reaction ). For different signes in eq.(13) the function 
has different properties with respect to crossing symmetry. For plus ( positive signature) 
the function is symmetric while for minus ( negative signature ) it is antisymmetric. 



2.2 s - channel unitarity 

To satisfy the s - channel unitarity we have to assume that trajectory a(t) < 1 in the 
scattering kinematic region ( t < ). This is why the exchange of the Reggeons can solve 
our first puzzle. 



2.3 Resonances. 

Let us consider the same function but in the resonance kinematic region at t > 4m 2 . Here 
a(t) is a complex function. If t — > t a(t ) = j = 2k where k = 1,2,3, .. the Reggeon 
exchange has a form for the positive signature: 

s 2k 

A R (s, t) 9l g 2 ■ — - _ - - Uma{tQ) (14) 

Since in this kinematic region the amplitude A R describes the reaction a + a — > b + b 

s = p 2 sin9 

where p = the amplitude has a form 

A ^ = 9i92 _ p 2j sin j 6 

r - gi92 ■ a ,( to ^ t _ to } + u ma (t ) ~ a '(t ) ' t - 1 + iT [ ' 

where the resonance width T = /m "fto) . 

a' (to) 

Therefore the Reggeon gives the Breit - Wigner amplitude of the resonance contribu- 
tion at t > 4m 2 . It is easy to show that the Reggeon exchange with the negative signature 
describes the contribution of a resonances with odd spin j — 2k + 1. 
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2.4 Trajectories. 



We can rephrase the previous observation in different words saying that the Reggeon 
describes the family of resonances that lies on the same trajectory a(t). It gives us a 
new approach to the classification of the resonances, which is quite different from usual 
SU 3 classification. Fig. 3 shows the bosonic resonances classified according the Reggeon 
trajectories. The surprising experimental fact is that all trajectories are the straight lines 

a(t) = a(0) + at (16) 

with the same value of the slope a' ~ lGeV~ 2 . 

We would like to draw your attention to the fact that this simple linear form comes 
from two experimental facts: 1) the width of resonances are much smaller than their mass 
(Tj <C MrJ and 2) the slope of the trajectories which is responsible for the shrinkage 
of the diffraction peak turns out to be the same from the experiments in the scattering 
kinematic region. 

The SUs means in terms of trajectories that 

a p (0) = a w (0) = 0^(0) 

in addition to the same value of the slope. The simple estimates show that the value of 
the intercept a(0) ~ 0.5 and therefore the exchange of the Reggeons give the cross section 
falling down as function of the energy without any violation of the Froissart theorem. 



2.5 Definite phase. 

The Reggeon amplitude of eq. (p~3|) can be rewritten in the form: 

A R = gig2V±s a(t) (17) 

where rj is the signature factor 

7ia(t) 

V+ = ct 9^~ + i 

Trait) 
V- = tg— % 

It means that the exchange of the Reggeon brings very definite phase of the scattering 
amplitude. This fact is very important especially for the description of the interaction 
with a polarized target. 
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2.6 Factorization. 



The amplitude of eq. ( |l3l) has a simple factorized form in which all dependences on the 
particular properties of colliding hadrons are concentrated in the vertex functions g\ and 
g 2 . To make it clear let us rewrite this factorization property in an explicit way: 



Ar = gi{m X: Afi, q 2 ) g2(m 2 , M 2 , q 2 ) • r]± ■ 



s 



For example this form of the Reggeon exchange means that if we try to describe the 
electron deep inelastic scattering with the target through the Reggeon exchange, only the 
vertex function should depend on the value of the virtuality of photon (Q 2 ) while the 
energy dependance does not depends on Q 2 . 



Problem 5: Show that the Reggeon exchange has the following form in bf. 

i i 

A R (s,b t ) = gi (p)g 2 (p)s«°).— j— ir — — -e ^^'^ 

An (Rf + Rein + a' Ins 

if we assume the simple exponential parameterization for the vertices: 

g 2 (q 2 ) = g 2 (0) e- R ^ 



2.7 Shrinkage of the diffraction peak. 

Using the linear trajectory for the Reggeons it is easy to see that the elastic cross section 
due to the exchange of the Reggeon can be written in the form: 



da, 



el 



dt 



9 2 M) ■ <li(<lt) ■ s 2(a(0) - 1} • e" 2Q 'W ln ^? (19) 



The last exponent reflects the phenomena that is called as the shrinkage of the diffraction 
peak. Indeed, at very high energy the elastic cross section is concentrated at values of 
q 2 t < ^- s . It means that the diffraction peak becomes narrower at higher energies. 
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3 Analyticity + Reggeons. 



Now we can come back to the main idea of the approach and try to construct the amplitude 
from the analytic properties and the Reggeon asymptotic at high energy. Veneziano 
[|IJ suggested the amplitude that satisfies the following constraint: it is the sum of all 
resonances in s - channel with zero width and simultaneously the same amplitude is the 
sum of the exchanges of all possible Reggeons. Taking the simplest case of the scattering 
of scalar particle, the Veneziano amplitude looks as follows: 

A = g[V(s,t) + V(u,t) + V(s,u)} (20) 

where 

v(s , t) = r ''-'^"7;;» (21) 

r(l — a(t) — a{s)) 

One can see that the Veneziano amplitude has resonances at a(s) — n + 1 where n = 
1, 2, 3..., since T(z) — > — \— at z — >• — n. At the same time 

' ' ' V / z+n 

A r(l-a(t))[(-«( S )) Q W + {-a{u)) a ®\ 

and therefore reproduces the Reggeon exchange at high energies. 

This simple model was the triumph of the whole approach showing us how we can 
construct the theory using analyticity and asymptotic. The idea was to use the Veneziano 
model as the first approximation or in other word as a new Born term in the theory and 
to try to build the new theory starting with the new Born Approximation. The coupling 
constant g turns out to be dimensionless and smaller than unity. This fact certainly also 
encouraged the theoreticians in 70 's to try this new approach. 



4 The second and the third puzzles: the Pomeron?! 

The experiment shows that: 

1. There is no particles ( resonances) on the Reggeon trajectory with the value of 
the intercept which is close to unity ( a(0) —*■!)■ As it has been mentioned the typical 
highest value of the intercept is a(0) ~ 0.5 which generates the cross section of the order 
of a tot oc s a ^~ x oc s~2. 

2. The total cross section is approximately constant at high energy. It means that we 
have to assume in the framework of the Reggeon approach that there is the Reggeon with 
the intercept close to 1. 
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The above two statement we call the second and the third puzzles that have to be 
solved in theory. The strategy of the approach was to assume that the Reggeon with the 
intercept close to 1 exists and tried to understand how and why it would be different from 
other Reggeons, in particular, why there is no particle on this trajectory. Now, let me 
introduce for the first time the word Pomeron. The first definition of the Pomeron: 

The Pomeron is the Reggeon with a(0) -1 e A < 1 

The name of Pomeron was introduced after russian physicist Pomeranchuk who did a 
lot to understand this miracle object. By the way, the general name of the Reggeon was 
given after the italian physicist Regge, who gave a beautiful theoretical arguments why 
such objects can exist in the quantum mechanics and the field theory 0. 

Let me summarize what we know about the Pomeron from the experiment ( see for 
example ref. ||): 

1. A ~ 0.08 

2. a'(0) ~ 0.25GeV~ 2 

Donnachie and Landshoff gave an elegant description of almost all existing experimen- 
tal data using the hypothesis of the Pomeron with the above parameters of its trajectory. 
However, we have to find the theoretical approach how to describe the Pomeron and why 
it is so different from other Reggeons. 

5 The Pomeron structure in the parton model. 
5.1 The Pomeron in the Veneziano model. 

As has been mentioned the Pomeron does not appear in the new Born term of our ap- 
proach. Therefore the first natural idea was to make an attempt to calculate the next to 
Born approximation in the Veneziano model to see can the Pomeron appear in it. The 
basic equation that we want to use is graphically pictured in Fig. 5, which is nothing more 
than the optical theorem. However we have to know the Born approximation for the 
amplitude of production of n particle. 

However to understand the main properties and problems which can arise in this 
approach let us calculate the contribution in equation of Fig. 5 of the first two particle 
state. This contribution is equal to: 

a tot = s «0)-i) J d 2 q' t T 2 {l-a{q' t 2 ))-e- 2a '^^ s (22) 
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since 

r(l - a{qf )) oc (Mlf^?) 

one can see that the essential values of q' in the integral is rather big, of the order of 
qf ~ s. It means that we have to believe in the Veneziano amplitude at the large value 
of momentum transfer. Thus the lesson that we have learned from this exercise is the 
following: 

To understand the Pomeron structure we have to understand better the structure of 
the scattering amplitude at large values of the momentum transfer or in other words we 
should know the interaction at small distances. 



5.2 The parton model. 

Feynman || and Gribov H suggested the simple model for the scattering amplitude at 
small distances that manifests itself in the deep inelastic scattering [§],so called the parton 
model. Let us assume that in the production amplitude the mean transverse momentum 
of the secondary particles does not depend on the energy ( ku = Const(s) ( see Fig. 6). If 
it is so, the main contribution to the equation of Fig.5 comes from very specific region of 
integration. Indeed, the total cross section can be written as follows: 

/dx ■ 
|M n 2 (x 4 p,^)|IL^ d 2 k tl (23) 
Xx 

where Xj is the fraction of energy that carried by the z-th particle. Let's call all secondary 
particle partons. 

It is clear that the biggest contribution in the above equation comes from the region 
of integration with strong ordering in Xi for all produced partons,namely 

2 

Tfl 

x 1 > x 2 > > Xi > x i+1 > > x n = — (24) 

s 

Integrating in this kinematic region we can put all Xi = in the amplitude M n . Finally, 

„„ = s„ /nA|A«M f\ ^ IT — - IT — (25) 

s fc s 6 s lb 

= S n |nA|M n 2 (^)| • ^ \n n s 
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This equation shows one very general property of the high energy interactions, namely 
the longitudinal coordinates ( Xj) and the transverse ones (k ti ) turns out to be sepa- 
rated and should be treated differently. In some sense the above equation reduced the 
problem of high energy behaviour of the total cross section to the calculation of the am- 
plitude M n which depend only on transverse coordinates. Assuming, for example, that 
/ Yld 2 kti\M^(k ti )\ 2 oc -^g n we can derive from the eq.(25) 

otot = \T, n 9 -\n n s = — 2 -s 9 (26) 
which looks just as Pomeron - like behaviour. 



Problem 6: Show that in g<p 3 - theory where is the scalar particle with mass m the 
total cross section is equal to 

g A 

&tot = — ir A — s47,m ' 2 

5.3 Random walk in b t . 

The simple parton picture reproduces also the shrinkage of the diffraction peak. Indeed, 
due to the uncertainty principle 

Ab ti hi ~ 1 (27) 

or in different form 

Ab u - ' 



<h> 

Therefore after each emission the position of the parton will be shifted on the value of 
Abt which is the same in average. After n emission we have the picture given in Fig. 7, 
namely the total shift in b t is equal to 

bl = r 2 ' n ( 28 ) 
tn < k t > 2 

which is typical answer for the random walk in two dimensions ( see Fig. 7). The value 
of the average number of emission n can be estimated from the expression for the total 
cross section, because 

i g n 



m? n\ m? nl 

and the value of < n >~ g Ins. If we substitute this value of < n > in the eq.(O) We get 
the radius of interaction 

R 2 = < b 2 n > — ^-^r - Ins = a' Ins 

< k t > 2 
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6 Reggeon Calculus. 



6.1 The main idea of the approach. 

The Reggeon Calculus was the first attempt to build the effective theory at high energy 
with the goal to define the asymptotic behaviour of the scattering amplitude. The brick 
from which we wanted to do this was the Pomeron. However, it turns out that the simple 
hypothesis that the Pomeron gives you the asymptotic behaviour of the amplitude is not 
correct. We have to take into account the interaction of the Pomerons. To illustrate this 
fact let us consider so called the triple Pomeron interaction which related to the process 
of the diffraction dissociation ( see Fig. 8 ). 

The cross section of the diffraction dissociation ( single diffraction SD) can be written 
in the form: 

where <jq is the total cross section at s = so due to exchange of one Pomeron and R\ = 
2Bsd, where Bsd denotes the slope of the SD cross section. It should be stressed that at 
very high energy when a' P > R%, R\ — > 4a' P In (s/M 2 ) 

Integrating over M 2 the SD cross section one can see that the total cross section of 

SD 

&SD OC ( — ) > a tot OC ( — ) . 
So s 

It means that we have to consider the interaction of the Pomerons to get the correct 
asymptotic at high energy. 



Problem 7: Show that energy s' in Fig. 8 is proportional to j^. 



Problem 8: Show that the total cross section of SD is larger that one Pomeron exchange 
even at A = 0. 
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6.2 Reggeon Diagram Technique. 

It turns out that the simplest way to deal with the Pomeron interaction is to use the 
Mellin transform of the amplitude |7j 

A(u,t) = / s'^ds ImA(s,t) (30) 
J o 



Problem 9: Show that the Mellin transform of the single Pomeron exchange gives 



u - A - a'pqf 
where t = —qf. 



Using the above Pomeron propagator it is easy to write any diagram for Pomeron 
interaction. For example, for the triple Pomeron interaction the diagram of Fig. 9 can be 
written in the form: 

1 1 

AsP = 9 7 —"I JTa TTT M 2 I 721 N ( 31 ) 

oj — A — a P qf to — 2A - oc P [ {q — k)i + kf \ 



6.3 The effective Lagrangian. 

The above examples show that the new effective theory can be build as the theory of the 
interaction of the effective degrees of freedom (Pomerons) in (1 + 2) dimensions. However, 
it is obvious that we cannot built the closed theory without specification what kind of 
interaction we want to take into account. On this pure phenomeno logical stage we have 
no selection rules and at first sight have to include any possible interactions. The idea was 
to start from triple Pomeron interaction, to solve the problem and to check back whether 
the more complicated interactions would be essential. 

For triple Pomeron interaction the effective Lagrangian looks very attractive: 

L = -\[^ + % + ^<f>] + - ^7[0 + + </#V (32) 
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where t = Ins and the propagator of the Pomeron is defined as the Green function of the 
Lagrangian without interaction (7 = 0). We will discuss a bit later the sign minus in 
front of the last term. 

The attempts to solve the effective theory were very interesting and we learned a lot, 
but we will not discuss them because we will show later that the whole approach is deadly 
sick. 

6.4 The AGK cutting rules. 

The AGK cutting rules || establish the generalization of the optical theorem on the case 
of the multi Pomeron interaction and in the simplest example of two Pomeron exchange 
they are pictured in Fig. 11. The meaning of the AGK cutting rules is very simple. Each 
Pomeron interaction contribution to the total cross section really corresponds to the cross 
section of the different processes. In the case of the two Pomeron exchange these processes 
are (i) the diffraction dissociation ( ) , (ii) the production of the secondary hadrons 
with the same multiplicity as for one Pomeron exchange ( o^ 1 -* ) and (iii) the production 
of the secondary hadrons with multiplicity in two times larger than for one Pomeron 
exchange ( ). 

The AGK cutting rules claim: 

a (o) + a (i) + a m = l + _ 4 + 2 (33) 

Two important consequences follow from the AGK cutting rules: 

1. The total cross section of the diffraction dissociation (DD) is equal to the contri- 
bution of two Pomeron exchange to the total cross section ( o^P ) with opposite sign: 

a (DD) = _ a (2P) 

2. The two Pomeron exchange does not contribute to the total inclusive cross section 
in central kinematical region. Indeed only two processes cW and are the sources of 
the produced particle in the central region. It means that the total inclusive cross section 
due to two Pomeron exchange is equal to 

Vine = CT (1) + 2 CT (2) = 

Factor 2 comes from the fact that the particle can be produced from two different parton 
showers (see Fig. 11). 
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Now let me give you a brief proof of the AGK cutting rules for the case of hadron 
- deuteron interaction. For simplicity let us assume that G in {b t ) — k — Const{b t ) for 
b t < R N and G in = for b t > R N , so the inelastic cross section for hadron - nucleon 
interaction is equal to 

cr% el = kS 

where S is the area of the nucleon ( S = tiR 2 n ). To calculate the elastic cross section we 
need to use the unitarity constraint in b t ( see eq.(7)) and the answer is 

<4 = 

Note that the flux of incoming particles after the first interaction becomes (1 — k) we can 
calculate the total inelastic interaction with the deutron, namely 

a™ 1 =kS + k(1-k)S = 2a™ el -n 2 S 

Using unitarity we have to calculate the elastic cross section for the interaction with 
deuteron which is equal to 

el ^ K \2q 

a D = ly) £> 

Therefore the total cross section for hadron - deutron interaction can be presented in the 
form: 



tot cy tot A 

r D = Za N — Act 



where 



2 ( inel\2 

Aa = Aa mel + Aa' ! -S [ ° " 



2 2ttR 2 

Everybody can recognize the usual Glauber formula for hadron - deutron interaction. 

As has been mentioned we have two sources of the inelastic cross section which are 
pictured in Fig. 13 for the case of the deutron. The cross section of the inelastic process 
with double multiplicity is easy to calculate, because it is equal to the probability of two 
inelastic interaction: 

A<xg } = k 2 -S 
To calculate;ate a$ we have to remember that 

Aa inel = a (2) + a (l) = _ kS 



Therefore = —2kS. Remembering that Acr^ — a% — 2a^ we get 
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n 





< n >n 


2< n > N 




K 2, Q 

2 D 


-2k 2 S 


K 2 S 



Therefore we get the AGK cutting rules for the hadron - deuteron interaction, since 
cr^* = 2<t^* corresponds to one Pomeron exchange and the correction to this simple for- 
mula just originated from the two Pomeron exchange in our Reggeon approach. However 
the above discussion, I hope, shows you that the AGK cutting rules have more general 
background than the Reggeon approach. For example, they hold in QCD providing the 
so called factorization theorem. 

Let me give here the general formula for the AGK cutting rules ||. If we have the 
contribution to the total cross section with the exchange of v Pomerons ( <Jt ot ), the cross 
sections of the process with the multiplicity of produced particles equal fi < n > which 
are generated by the same diagram (see Fig. 13) are equal[|: 

— Uo = -ir" • -n T7 • 2" (34) 

while for /i = the ratio is equal 

— = (-ir-[2^ 1 - 1] (35) 

°tot 

For the first 5 exchange you can find these factors in the following table. 
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- 3 
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- 12 
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- 32 


48 


- 32 
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- 15 


80 


- 160 


160 
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6.5 Different processes in the Reggeon Approach. 

The AGK cutting rules together with Mueller theorem establish the relationship be- 
tween the contribution of many Pomeron exchanges and different exclusive and inclusive 
processes in the Reggeon Approach. Here we want to write down several examples of 
different processes that can be treated on the same footing in the reggeon approach. 

* < n > is the average multiplicity in one Pomeron exchange 
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6.5.1 Total cross section. 

As has been mentioned in one Pomeron exchange approximation the total cross section 
is given by the following expression ( see Fig. 15 []): 

dot = 47^(0) s 2 (0)(-) A a (-) A (36) 

So So 

We would like to remind you that multi Pomeron exchange is very essential in the total 
cross section but we postpone the discussion of their contribution to the second part of 
our lectures. 



6.5.2 Elastic cross section. 



Problem 10: Show that for one Pomeron exchange the total elastic cross section is equal 
to 

O 'tot 



where 

B el = 2R 2 01 + 2R 2 02 + 2a' P lns/s 
in the exponential parameterization for the vertices g(t) ( see Problem 5). 



6.5.3 Single diffraction dissociation 

The cross section of the single diffraction (see Fig. 15) has the following form: 



• 

Problem 11: Show that in exponential parameterization for vertices (see Problem 5) 

R\^) = 2Rl 1 + rl + 4a' P Hs/M 2 ) 

tin Fig. 15 we use a little bit different notation, namely (3(t) instead of g (t ) and g ( t) for G$p{t). 
We did this because we had nice prepared pictures in such notations. Hope that this fact will not lead 
to misunderstanding. 
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where r 2 . is the radius of the triple Pomeron vertex which experimentally is very small 
(r 2 < lGeV- 2 ). 



One can see that we include in the expression for the cross section of the single diffrac- 
tion also the contribution of the secondary Reggeon trajectory which describes the be- 
haviour of the cross section at rather small values of produced mass (M 2 ). 

6.5.4 Double diffraction dissociation. 



Problem 12: Show that the cross section of the double diffraction dissociation ( DD 
) process (see Fig. 15) is equal to: 

MlMld 2 o DD _ a 2 ssp a2A M| A M[ A 



in the region of large values of produced masses ( Mi and M 2 ). 



Problem 13: Show that Rq in the above equation is equal 

Rl = Kl + Aa' P \n{ss G /MlMl) 
in exponential parameterization of the vertices ( see Problem 5 ). 



Comparing the cross section for double and single diffraction as well as for elastic cross 
section one can get the following factorization relation: 

,^9,^9,9 Mfdasp M%da SD p2 / s \r>2/ s \ 

dMldM 2 : a el ' R 2 o(j^)R 2 el (s) [ ] 

where 

^ = 2^ + 2^ + 4^^(5/50) 

The structure of the event which corresponds to double diffraction is pictured in Fig. 15 
in lego - plot. In the region of rapidity Ay we have no produced particles. 
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6.5.5 Central Diffraction. 



In this process we produced the bunch of secondary particle in the central rapidity region, 
while there are no secondary particles in other regions in rapidity. 

The first diagram that describes this process gives the answer: 

Kq 

The above examples demonstrate how we can calculate the exclusive processes using 
several phenomenological inputs in the Reggeon Approach, however all these processes 
are crucially affected by multi Pomeron contributions and this subject we are going to 
discuss later. 



6.5.6 Inclusive cross section. 

As has been discussed the inclusive production does not depend on multi Pomeron ex- 
changes. Using Mueller technique (see Fig. 16) which is a generalization of the optical 
theorem for more complicated than total cross sections cases the inclusive cross section 
can be written in a very simple form: 

-:- = a-a tot (41) 

dy c 

where a is the new vertex for the emission of the particle c. The inclusive reaction that 
is under consideration is 

a + b — > c(y) + anything 



6.5.7 Two particle rapidity correlation. 



The Reggeon approach we can use for the estimates of the two particle rapidity correlation 
function, which is defined as: 



d?(r(yi,y2) 
atotdyidm 
Irr(yi) da(y 2 
dyi dy 2 



r = . -i (42) 
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where , f is the double inclusive cross section of the reaction: 

dy\dyi 

a + b — > l(yi) + 2(y 2 ) + anything 



Problem 14: Show using the AGK 
equal: 

i?(Ay = |yi-y 2 |) = 
All notations are clear from Fig. 16. 

It should be stressed that the Reggeon approach gives the estimates for the correlation 
length ( Lcor) and the the strength of the long range correlations ( note, that the second 
term does not depends on rapitidies of produced hadrons). Namely, we can rewrite the 
above formula in the form: 

Ay 

R = SR-e~— + LR 

and from our Reggeon formula we have L cor ~ 2. To evaluate the second term we have 
to develop some model for including the multi Pomeron exchanges in the calculation of 
the total cross section. We will do this in the second part of the lectures. 



cutting rules that the correlation function R is 



,2 J2P) 



PR . e (l-a«(0))A„ + 2 ^o* (43) 

pp (Hot 



6.6 Hopes. 

The main idea of the Reggeon Calculus was to build the effective theory for the strong 
interaction at high energy starting from the effective Lagrangian and the AGK cutting 
rules. The last one played a very crucial role giving us the possibility to look insight of 
the inelastic interaction using only general ideas on the structure of the Pomeron. 

Of course one defect of the approach has been seen from the beginning, namely the 
absence of a theoretical idea how to select the interaction between Pomerons. However 
we could hope that the future theory will provide us such selection rules and we will be 
able to adjust the developed formalism. 
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7 The death of the Reggeon Approach. 



In 1974 - 1975 in ref. |L0[ was shown that all our hopes were in vain. T.T. Wu and 
McCoy and S.Matinian and A. Sedrokian proved that in the general class of the theories ( 
all theoretical model that existed at that time) our hope that the Pomeron and Pomeron 
interaction can describe the high energy asymptotic is not correct. 

Let me illustrate this point using the parton model. In this model the Pomeron is the 
sum of "ladder" diagrams while the two Pomeron exchange corresponds to the diagram 
of Fig. 17 where the partons are produced and absorbed by the same "ladders". These 
diagrams have the probabilistic interpretation and look very natural in the parton model. 
However it turns out that the diagrams where the secondary partons are produced by one 
"ladder" and are absorbed by the second one give bigger contribution. These diagrams 
have not been included neither in one Pomeron exchange nor in two Pomeron exchange. 

It means that our hope to reduce the whole theory to Pomeron interaction can not be 
right or at least we have to build more complicated picture for the Pomeron structure, 
which is impossible to do without detailed microscopic theory. 

It means also that the AGK cutting rules we cannot apply for the process of the 
Pomeron interaction. For example, if we want to calculate the inclusive production of the 
hadron in the single diffraction process we cannot restrict ourselves by calculation of only 
two diagrams of Fig. 18, but we should add also the third diagram. 



8 Lessons for future. 

Inspite of the fact that we failed to construct the effective theory for high energy inter- 
action starting from the new effective degree of freedom -Pomeron we certainly learned a 
lot about properties of high energy interaction. In this section I want to summarize those 
lessons that we have learned and which could be useful for future. 

1. The longitudinal and transverse degrees of freedoms look differently at high energy 
and can be treated separately and in different ways. 

2. The effective theory at high energy can be reduced to two dimensional theory and 
we need to find the energy spectrum of this two dimensional theory to specify the high 
energy asymptotic. 

3. The major problem of any affective theory at high energy is to find correct degree 
of freedom or in other words the correct effective particle which is responsible for high 
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energy asymptotic. The concrete realization, namely Reggeon Calculus failed, because the 
Pomeron, effective particle of that time, turns out to be not the correct one to construct 
theory at high energy ( see the previous section). 

4. The only way to invent correct effective particle is to develop the consistent approach 
starting from microscopic theory but not from phenomenology. 

At that time we had no microscopic theory and in the best tradition of high energy 
physics the experts left the field. Now situation is quite different, we have good microscopic 
theory ( QCD) and certainly we have a lot of problems in QCD which have to be solved. 
High energy asymptotic is only one of many. I think it is time to ask yourselves why we 
spend our time and brain trying nevertheless to find the high energy asymptotic in D1S. 
Let me give the first answer to this question, namely it seems extremely interesting to go 
back to old problems to see how far away was our guess from the systematic approach. 
It is rather private but strong motivation for many experts including me. Later we will 
discuss this problem in more details. 

We want to finish our review of the past in high energy physics with a statement, 
that inspite of all theoretical inconsistency, all intristic contradictions of the Reggeon 
Calculus, this old fashioned approach is still the main source of our terminology in high 
energy theory and still the only phenomenological approach that we have in hand even 
now when we discuss so called "soft" processes. It is why the knowledge of this approach 
is still the alpha and omega of high energy strong interaction theory. We hope that we 
gave you just right amount of information to move to more theoretical approach based 
on QCD. 
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TODAY: 



POMERON WITHOUT MAGIC 
IN PERTURB ATIVE QCD 

9 Basics of perturbative QCD (pQCD) 

9.1 The basics of pQCD for high energy scattering. 

In this subsection we are going to discuss the main ideas of perturbative QCD, which is our 
microscopic theory of the strong interaction. The goal of our approach is to understand 
what is the Pomeron in QCD, but we start from the several things that everybody knows 
about QCD. 

1. If the typical distances (r) in our process are small we have a natural small param- 
eter in QCD, namely the strong coupling constant 

/ \ 4?r «s(/i 2 ) ^ i 

<*<'> = = 1 + « 1 

where b — 11 — rif is the number of quarks, A is the QCD confinement scale and /i 
is the renormalisation scale. 

It means that we can start to calculate the amplitude of our process using the expansion 
with respect to small a$. 

2. The Born approximation of pQCD gives the cross section which is constant at high 
energy 

a t — > | s ^oo Const 
We will discuss this property in details in the next section. 

3. The more complicated diagrams lead to increase of the total cross section with 
energy as we will discuss later 



where A = Ca$. 
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All the above properties show us that we have a good chance to get in QCD the high 
energy asymptotic which looks like the one Pomeron exchange ( the word Pomeron I use 
here in the traditional sense as the Reggeon with intercept close to 1 ) and we are able 
to examine how good or bad is our traditional Reggeon - like approach. However, we are 
going to start with the brief review what we have learned about pQCD to give you an 
impression what kind of theory we have to attack the Pomeron structure. 

9.2 The map of QCD. 

Describing the basic features of QCD, we will start to answer in this subsection the 
following questions: (i) What we have learned about QCD; (ii) What problems are still 
unsolved in QCD?; (iii) Can we relate the kinematical region in which we have to face 
these unsolved problems to collision processes? 

Fig. 19 displays the 'map of QCD'. Shown are three separate regions, distinguished 
by the size of the variables Q 2 and x, which will be defined momentarily. Each region 
corresponds to quite different physics and enjoys a different level of understanding. Before 
we discuss these three regions in more detail let us introduce the necessary notational 
conventions and definitions. 

Deeply Inelastic Scattering ( DIS ). 

In Fig. 19 r is the distance that can be resolved by the scattering process under consid- 
eration. From the uncertainty principle this distance is of the order of ^ where Q is the 
typical large transverse momentum in our experiment. E.g. in deeply inelastic electron 
scattering Q is the transverse momentum of the recoiled electron. We recall that deeply 
inelastic scattering is the reaction: 

e + p — > e + anything. 

Thus, this reaction acts as a powerful microscope which is able to resolve the constituents 
of a hadron (quarks, antiquarks and gluons, collectively called partons) with a transverse 
size of the order of ^. The second kinematical variable that we can introduce for such 
constituents is the fraction of energy (x) that a parton carries with respect to the parent 
hadron. 

It should be stressed that the energy in DIS is equal to s = It means that we 
are interested in low x DIS if we are going to discuss the high energy asymptotic. If 
N(x, Q 2 )dx is the number of gluons in a small x-interval centered around the value x at 
scale Q 2 . We will show later that the number of gluons increases in the region of small x. 
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All physics in this kinematical region is strongly connected to this fact. It is the reason 
we concentrate on the discussion of the gluon density here. 

The gluon structure function xG(x, Q 2 ) tells us what the gluon density at a definite 
value of In - is, 

2 i dN(x,Q 2 ) i 
xG(x,Q) = \—^-\. 

X 

Let us for convenience introduce the transverse density of gluons p, 

p(x,Q 2 ) = X -^2^1 (44 ) 

where R is the radius of the hadron. 

We plot this on the vertical axis in Fig. 19. In this figure we can see three different 
regions: 

1. The region of small parton density at small distances (low density (pQCD) region). 

This is the region where we can apply the powerful methods of perturbative QCD since 
the value of running coupling constant a s (l/r 2 ) is small (a s (l/r 2 ) <C 1). During two 
decades remarkable theoretical progress has been achieved here (GLAP evolution equation 
gluon bremshtrahlung for jet decay |12 , factorization theorem (J. Collins, D.Soper 



and G.Sterman (1983) and the main property of "hard" processes has been experi- 
mentally confirmed at LEP and at the Tevatron. 

2. The region of large distances (npQCD region). 

Here we have to deal with the confinement problems of QCD, since a s (l/r 2 ) ^> 1. In 
this kinematical region we need to use nonperturbative methods such as lattice calculation 



or QCG Sum Rules |T5|. The progress here is remarkable but all developed methods 



cannot yet be applied to scattering processes. 

3 The region of small distances but high parton density (hdQCD region). 

Here we have a unique situation in which the coupling constant a s is still small but 
the density is so large that we cannot use the usual methods of perturbation theory. 

In essence the theoretical problem here is also a nonperturbative one but the origin of 
the nonperturbative effects here is quite different from that in the previous region. Here 
we face the situation where we have to develop new methods that let deal with a dense 
relativistic system of gluons in a nonequilibrium state. Unfortunately we are only at the 
beginning of this road. 
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Fortunately, we can control theoretically this dense system of partons in some transi- 
tion region on the border of the pQCD and hdQCD regions and here we can study this 
remarkable system of partons in great detail. Thus the right strategy is to approach this 
interesting kinematic region from the low density pQCD region. 

In my opinion the interest to high energy scattering was revived only because moving 
from pQCD region to hdQCD region we can open new window to study a nonperturbative 
nature of QCD. 

9.3 What fundamental problems could be solved? 

This is why we want to list here the fundamental problems that we hope to solve pene- 
trating high density QCD region. We hope: 

1. to specify the kinematical region in which we can trust pQCD (GLAP evolution 
equation, gluon bremstrahlung, factorization theorem ...); 

2. to find new collective phenomena for nonabelian theories such as QCD ; 

3. to find the analytic solution of hd QCD which is nonperturbative but looks simpler 
than np QCD since a s -C 1 here; 

4. to develop methods with which build an effective theory for hd QCD. 

We need the effective theory because we can use the Lagrangian of such a theory for 
exact calculation on lattice, for example, or even we can try to solve problem analytically. 

9.4 How to penetrate the high density QCD region. 

Access to this interesting kinematical region is actually easily achieved in our scattering 
processes. We know at least three ways to prepare a large density system of partons. 

1. The first is given by nature, which supplies us with large and heavy nuclei. In ion- 
ion collisions we can already reach a very high density of partons at not so high energies, 
because the partons from different nucleons in a nucleus are freed. 

2. The second relates to hard processes in hadron-hadron collisions or in deep inelastic 
scattering. These also give us access to a high density of partons because we expect a sub- 
stantial increase in the parton density in the region of small Bjorken x. The experimental 
data from HERA show the significant increase of the deep inelastic structure function: 

F 2 (Q 2 ,x B ) oc (— ) a33 atQ 2 ~ WGeV 2 . 
x B 
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3. The third is to measure the event with sufficiently large multiplicity of produced 
particles, larger than the multiplicity in the typical inelastic (bias) event. 

Thus we can formulate the ideal experiment for search of high density partonic system: 

" The deeply inelastic scattering with nucleus at low x with special selection of events 
with large multiplicity of produced particles. " 

9.5 Method of pQCD( two diabols that we are fighting with). 

Accordingly to our strategy we start with perturbative QCD in that region of Fig. 19 
where the density and coupling QCD constant (a s ) are small. Each observable (i.e gluon 
structure function) could be written in pQCD as following series: 

xG(x, Q 2 ) = S„ =0 C n (a s ) n ■ (L n + a n . x L n ~ x ...o ) , . (45) 

We have two big problems with this perturbative series which are two our biggest enemies: 

1. The natural small parameter a s is compensated by large log (L). The value of L 
depends on the process and kinematic region. For example in deeply inelastic scattering 
(DIS): 

A = log Q 2 at Q 2 ->Q 2 but x ~ 1 

A = \og(l/x) at Q 2 ~ Ql andx — > 

A = logQ 2 • \og(l/x) at Q 2 > Ql andx -> 

A = log(l — x) at Q 2 ~ Ql and x — > 1 

Of course it is not the full list of scales. The only that I would like to emphasize that L 
depends on the kinematic region. Thus to calculate xG(x, Q 2 ) one cannot calculate only 
the Born Approximation but has to calculate the huge number of Feynman diagrams. 

2. 

Cn — > n\ at n ^> 1 

It means that we are dealing with asymptotic series and we do not know the general 
rules what to do with such series. There is only one rule, namely to find the analytic 
function which has the same perturbative series. Sometimes but very rarely we can find 
such analytic function. In this case this is the exact solution of our problem. Mainly we 
develop some general approach based on Leading Log Approximation (LLA). The idea is 
simple. Let us find the analytic function that sums the series: 

xG(x,Q 2 ) LLA = Z n=0 C n (a s -L) n . (46) 
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Ussualy we can write the equation for function xG(x, Q 2 )lla- The most famous one, 
the GLAP evolution equation [11], sums eq.(2) if L— logQ 2 . The BFKL equation 



gives the answer for eq.(2) in the case when A=log(l/x). Using the solution of the LLA 
equation we built the ratio: 

R(x, Q 2 ) n^^? = S -° r " = E -=oCn • i^' 1 + *n-iL n - 2 + ...a ) . (47) 

This ratio is also asymptotic series but what we are doing we calculate this series term 
by term. Our hope is that the value of the next term will be smaller then the previous 
one (7-^ ^ 1) f° r sufficiently large n. However we know that at some value of n = N 
- J ^ L - ~ 1. The only that we can say about such situation in general that our calculation 
has intristic theoretical accuracy and the result of calculation should be presented in the 
form 

R(x,Q 2 ) = K=o~ X r n ±r N . (48) 

How big is the value of N depends mostly on how well we chose the LLA and how well 
we established the value of scale L in the process of interest. 



9.6 Present Theoretical Status = Regeneration of Reggeon Cal- 
culus. 

What we are doing now approaching the high density QCD domain is really the regener- 
ation of the old idea of Reggeon Calculus, namely we reduce the complicated problem of 
quark and gluon interaction in the dense parton system to the interaction of our "building 
bricks", so called Pomerons ( see Fig. 20, where the structure of our approach is shown). 
It looks like in old, good time of Reggeon Dominance but we have two new and very 
important ingredients: 

1. The QCD Pomeron is not an invention but naturally appears in perturbative QCD 
in the leading log(l/x) approximation (LL (x)A) to the scattering amplitude at high 
energy ( the Balitski - Fadin -Kuraev - Lipatov (BFKL ) equation [|16|] ). It means that in 



a restricted kinematical region the BFKL Pomeron describes the high energy interaction 
within a certain guaranteed theoretical accuracy. This fact makes it unavoidable that one 
should build an effective theory starting with the BFKL Pomeron. 

2. The vertices of the Pomeron interactions are not phenomenological parameters but 
can be calculated in perturbative QCD. It turns out that only two of them are essential 
in the vicinity of the border between pQCD and hdQCD regions in Fig. 19, namely triple 
Pomeron vertex (7) and rescattering of two Pomerons (A) (see Fig. 20). 
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We are going to discuss in detail the structure of the BFKL Pomeron as well as the 
theory of the shadowing corrections. The process which we consider is the deeply inelastic 
scattering in which all high energy phenomena occur at small distances where we can trust 
the pQCD calculation. 



10 Pomeron in pQCD (the BFKL Pomeron) 

The BFKL equation was derived in so called Leading Log (1/x) Approximation, in which 
we would like to keep the contribution of the order of (a$log{l / x)) n and neglect all other 
contributions, even of the order of as log(Q 2 /Qq). So the set of parameters in LL(l/x)A 
is obvious: ^ 

a s log- ~ 1 ; 

x 

Q 2 

a s log-^ < 1 ; (49) 

«s < 1 ; 

Let us consider the simplest process: the quark - quark scattering at high energy at zero 
momentum transfer . All problems of infrared divergency in such a process is irrelevant 
since they are canceled in the scattering of two colorless hadrons . We will show this fact 
using the example of the first diagram (the Born Approximation (BA)). 



10.1 The Born Approximation 

In the Born Approximation the only diagram of Fig. 21(a) contributes to the imaginary 
part of the scattering amplitude (A). It is easy to understand that the result of calculation 
of this diagram gives [|: 

2Im{^( M = 0)} = / J^l M(2 2; a s | s,t = -*») |> = S §± f **S , 

' (50) 

where M(2 — > 2; as | s, t = —kf) denotes the amplitude in the lowest order of as for 
quark - quark scattering at transfer momentum t = —k 2 through one gluon exchange ( 
see Fig. 2 1(a) ), N is the number of colours and C 2 = One can recognize the Low - 

Nussinov mechanism [|T^] of high energy interaction in this simple example. 

■fLet me recall you that we are keeping the coupling constant fixed. 
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Let us consider the BA in more detail for the scattering of two mesons ( a and b ), 
assuming that each of them consists of two heavy quarks. Such a system has sufficiently 
small size, namely 

1 

r oc — 

M Q a s 

and can be treated in pQCD. 

Let me remind you why the size is small. The total energy of two quarks in Coulomb 
potential is equal: 

E . . J? 

M Q r 

If we remember the uncertainty principle 

1 

ApAr ~ 1 or p ~ - 

r 

we get for the system with E — > that size which we have mentioned before. 

Repeating the calculations we get the same answer as for quark - quark scattering, 
namely 

""9 ' '^ 6 7 (fc - §)? (* + f )? 

where ra a (ra&) is the number of quarks in the meson a(6), the factor | is the result of 
averaging over colours of the initial quarks and $(&;, q) is the amplitude of the emission 
of two gluons by a meson. 



2 r d 2 k t 

a" = I - • 1 • "„"/, y , fc _ g x 2 / fc + £N2 ' ®a(k,q) $b(k,q) (51) 



Problem 15: Show that the emission amplitude $ can be expressed in terms of meson 
formfactors if we can treat the heavy quarks as nonrelativistic particles (see Fig.22(c)) : 

G(q 2 ) a = fdpy a (p)V:(p-q); 
$ a (k,q) = Jdpy a (p)y* a (p-q) - Jdp^ a (p+^ + k)^:(p+^-k) = G a {q 2 )-G a {Ak 2 ). 



Finally, we have 

2 4 • n a n b / • [ G a (0) - G a (4k 2 ) ] [ G b (0) - G b (4k 2 ) ] (52) 



a = ao 
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Now let us ask the question what kinematic region of integration over k t gives the biggest 
contribution to the total cross section if two mesons have different sizes, let say r a ^> r b . 
We can study three regions of integration: 

1. A > r a , where A is the wave length of the gluon (A = ^ ). 

In this case k t r a < 1 as well as k t r b < 1 and 



4 



2 7,2 



G a (0) - G a (4k ) - -r 2 a k 



G b (0) - G b (4k 2 ) - \r\kl 



Q a-t 

4 

-J 

6 



The integral over k t gives you the answer 



2.r a > A > r b 

In this case k t r a > 1 as well as k t r b < 1 and 

G a (0) - G a (4k 2 ) ^ 1 

G b (0) - G b (4k 2 ) - Uk? 



6 



The integral over k t gives you the answer 



a BA = al^-4-n a n b .l.rl\n r 4 (53) 
9 3 r b 

3. r b > X 

It means that k t r a > 1 as well as k t r b > 1 and 

G o (0) - G a {Ak 2 ) ^ 1 
G b (0) - G b (4k 2 ) ^ 1 

The integral over k t gives you the same answer as in the first kinematic region, namely: 

a BA ~ r 2 b 

This example is very instructive, because our standard process - DIS which we can treat 
using pQCD and in DIS the typical size of our probe - virtual photon is the smallest 
one in the process and is of the order of 4. It means that the BA gives for the DIS the 
total cross section which is proportional ^ ln(Q 2 R 2 ) where R is the size of the target. It 
should be stressed as has been mentioned that the total cross section does not depend on 
energy in BA. 
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10.2 The next order of as approximation. 



In the next order we have to consider a larger number of the diagrams, but we can write 
down the answer in the following general form: 



/ 



2Im{A" BA (s,t = 0)} = (54) 
(2vr)V 4 )( Pl + p 2 - P ; - p' 2 - q)| M(2 - 3; g 3 ) | 2 n^-^|- + 



J (2tt) V 4 >( Pl + Pl - pi - p' 2 ) ■ 2 Re{M(2 - 2; a 2 ) . M *(2 -> 2; a s )}Il£? 



(2vr) 3 2£; 



where g is coupling constant of QCD ( as = f-); M(2 — > 3;g 3 ) is the amplitude for 
the production of the extra gluon in the Born Approximation, and is given by the set of 
Feynman diagram in Fig. 22 while M (2 — > 2, a|) is the amplitude of the elastic scattering 
in the next to leading Born Approximation ( see Fig. 23) at the momentum transfer kt- 



Two terms in eq. flM ) have different physical meaning: the first one describes the 



emission of the additional gluon in the final state of our reaction, while the second term 
is the virtual correction to the Born Approximation due to the emission of the additional 
gluon. It corresponds to the same two particle final state, and describes the fact that due 
to emission the probability to detect this final state becomes smaller ( we will see later 
that the sign of the second term is negative). 



In the both terms of eq. fl5"4"D we can integrate over p 3 ' as well as over the longitudinal 
component of pi (p'ix). Finally, we rewrite the phase space in the following way: 



/(27r)V 4 )( Pl + p 2 -pi-p 2 -qniz?- 



d 3 p'i 1 1 f 1 dx' 3 f d 2 p' lt d 2 p' 3t _ 



2tt) 3 2^ 4vr s h mm =^- x' 3 J (2tt 



/(2^>( P1 + Pl - pi - py • m ^« = i ■ / g« ; (55) 



null 
2| 



(27r) 3 2E'i s J (2tt) 

where x is the fraction of the longitudinal momentum carried by particle, the value of x, 
depends on the reaction. For example in the deeply inelastic scattering x min = xb = 

2 

In the case of the quark scattering x m i n = ^f- where m t is the transverse mass of produced 
quark. 



From eq. (|55| ) one can see the origin of the log(l/x m i n ) contribution: it stems from 
the phase space integration, if M(2 — > 3) does not go to zero at x 3 —>■ 0. To sum 
the diagrams of Fig. 22 in this limit we can use two tricks. The first one is for each t- 
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channel gluon one can rewrite the numerator of the gluon propagator at high energy in 
the following way: 



9^u = — — + — 56 

pi • p 2 s 

The second trick is based on the gauge invariance of the QCD. We look on the subset of 
the diagrams of Fig. 22 pictured in Fig. 24 as the amplitude of the interaction of gluon 
k with the quark p 2 ( see Fig. 24 ). Since all particles in amplitude M v except gluon k are 
on the mass shell, the gauge invariance leads to the relationship: 

k v M v = . (57) 

Using Sudakov variables []TS1 we can expand vector k as 

k v = a k piu + fikVin + k tf M 

and rewrite eq. (|57D in the form: 

( a kPlu + (3 k p 2fl + k tl , ) M v = . (58) 

We note that all particle inside M v have a large component of their momentum on p 2 . It 
means that we can neglect the projection of vector M v on pi, or in other words 

M v = M^p lu + M^p 2v + AfW 

and M^> . Thus we can conclude from eq. (]58f) that 

p lv M v = - . 59 



Using both tricks of eq. fl56|) and eq. (|59D one can easily see that only diagram of Fig. 
24(1) contributes in LL(log(l/x)A. Indeed, let us consider for example the diagram of 
Fig. 24(2), the dominator of the quark propagator (p 2 + k) 2 is equal to 



72 m /2 



(pa + kf = a k s + k\ = + k\ ~ . 

Since due to eq. (^) the polarization of gluon k is transverse we cannot compensate the 
smallness this diagram at x% — > 0. Using the same tricks with the upper parts of the 
diagrams of Fig. 22 we arrive at the conclusion that that the set of the diagrams of Fig. 
22 degenerates into one diagram of Fig. 25 with specific vertex for gluon emission: 

T a = ig f abc 7^ , (60) 

akPk'S 
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where is given by the usual Feynman rules for QCD. 



Substituting eq. ( pOj ) into the first term of eq. we get the contribution of the 
emission of one additional gluon to the next to Born Approximation in the form: 

ImK™ „} = . / o " dy'J^-^. K cm ,„Uk„ K) ■ , (61) 

where y' = log(l/xz), y = log(l/x m i n and the kernel K(kt,k' t ) is equal to 



1,2 U2 

h - Kf 



emissiorjj^ti ^t) / 1 / \2 (^) 



To calculate the virtual correction in the next to Born Approximation (M( 2 — > 2 ; a| ) 
) we have to estimate the contribution of the set of diagrams of Fig. 23. The log (1/x) 
contribution is hidden in the real part of the amplitude M{ 2 — > 2 ; a| ), and easiest way 
to extract this log is to use the dispersion relation: 



TX J s — s u — u 

(63) 

We calculate Im{M(2 — > 2; ag)} s and Im{M( 2 — > 2; o|)}m using unitarity ( see 
eq. ©): 

Im{M(2 -> 2; ajj \ t = -k 2 t )} s = (64) 
/ Re { M ( 2 ^ 2 - «s I s,* = -fc?) M(2^2; a s I s, t = -(fc* - ^) 2 )} . 



The difference between eq. (|100|) and eq. (p0|) is that the Born amplitude for one gluon 
exchange enters these two equations at different values of momentum transferred t. The 
explicit calculations give: 

Im{M(2 -> 2; al \ t = -k 2 t )} s = s ■ C s vr£(A; 2 ); (65) 

Im{M(2 ^ 2; a 2 s \t = -k 2 )} u = s ■ C u ttE(A; 2 ); 



7T^ 



2 J {k t -k' t fk? ' 



Using the dispersion relation of eq. (|63|) one can reconstruct the real part of the amplitude 
and the answer is 

Re{M( 2 ^ 2; a 2 |t = -A; 2 )} = s ( C u - C s ) ■ S(A; t 2 ) • logs (66) 



37 



The colour coefficients have a very famous relation between them ( see Fig. 26 ) which 
gives, for the difference of the colour coefficients in eq. (|66|) the same colour structure as 
for the diagram of Fig. 23 ( 3 ). Thus Re{M( 2 — > 2 ; a 2 , 1 1 = —k 2 )} has the same 
colour structure as one gluon exchange in the Born Approximation. This fact makes it 
possible to rewrite the second term in eq. fl5"4] ) as the correction to the gluon trajectory, 
e.g. instead of gluon with propagator p • s we can introduce the new propagator 

^ • s a °^ (67) 
k 



a s N r k 2 t d 2 k' t a s N r k 2 t d 2 k' t 

~^~'J (kt - k' t yk' 2 ~ ~w'J {(k t - k' t y + k' 2 ]k' 2 



The answer for the second term in eq. fl51| ) can be written in the form 



Im{ A™ ,} = » #^7 / ( « G (* 2 ) - 1 ) • 2 • ^ • (® 



We can get the full answer for the amplitude in the next to the Born Approximation (a| 
by summing eq. fl61~|) and eq. ( |68"D and it can be written in the form: 

lmint = 0)} = sA^L r dy , f . z» . K{klX) . *k m 



N 2 — 1 Jo y J k 2 t vr 2 v 11 k[ 2 
where y = log(l/x min and 

111 k 2 1 

K{kuK) ' W = - k' t ) 2 ' V 2 ~ (*t-*0 2 [(*t'- + ' ¥ ' (70) 

Using eqs. (|69|)-([70|) we can introduce function <fi(k 2 ) and rewrite the total cross section 
for quark - quark scattering in the form: 

a$C 2 f , ( , 2 \ dk 2 

a w = N 2 _ 1 -J ^(y^)--jg- • ( 71 ) 



For eq. (|69f) gives the equation: 



(2) (V, k 2 t ) = f V dy' [ d 2 k' t K(h, K) (V', K) • (72) 



n 2 



where 



= \^ 2 -^ 2 ) - 77 — rj^ 1 ,m 2 i r^r -^ 2 )- (73) 



and 

/AW 



fc 2 



2 



(74) 
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10.3 The main property of the BFKL equation. 



From the simplest calculation in «g order one can guess the BFKL equation [IB| which 
looks as follows: 



d(f)(y = log(l/x),kf) a s N 



dy 7i 



2 



K( k t , k' t ) 0( y, k' t 2 ) d% , (75) 



where kernel K(k t , k' t ) is defined by eq. (|73|). This equation sums the (aslog(l/x)) n contri- 
butions and has "ladder" - like structure ( see Fig. 27 ). However, such "ladder" diagrams 
are only an effective representation of the whole huge set of the Feynman diagrams, as 
explained in the simplest example of the previous subsection. The first part of the kernel 
K(kt, k' t ) describes the emission of new gluon, but with the vertex which differs from the 
vertex in the Feynman diagram, while the second one is related to the reggeization of all 
t-channel gluons in the " ladder" . 



The solution of the BFKL equation has been given in ref. |TB| and we we would like 
to recall some main properties of this solution. 

10.3.1 Eigenfunctions of the BFKL equation. 

The eigenfunction of the kernel K(kt,k ; t ) is </>/ = (kf)^~ x . Indeed after sufficiently long 
algebra we can see that 



d 2 k' t K(k t ,k' t ) MK 2 ) = x(f) MW ( 76 ) 



i 

7T 

where 

X (f) = 2*(1) - *(/) - (77) 

and 

*(/) - ^> , 

r(/) is the Euler gamma function. 

10.3.2 The general solution of the BFKL equation. 

From eq. (|76|) we can easily find the general solution of the BFKL equation using double 
Mellin transform: 
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where the contours of integration over oo and / are situated to the right of all singularities 
of 4>(f) and C(u>, f). For C(w, f) the equation reads 

uC(u,f) = S**L x (f)C{u,f). (79) 

71 

Finally, the general solution is 



^y,k?) = / dfe^^y + ^ r cf>(f) (80) 

where <p(f) should be calculated from the initial condition at y = yo and r = ln^ ( q$ 
is the value of virtuality from which we are able to apply perturbative QCD ) . 



10.3.3 Anomalous dimension from the BFKL equation. 



We can solve eq. ( [79] ) in a different way and find / = j(oo). 7(w) is the anomalous 
dimension in LL(log (l/x)A | and for 7(0;) we have the following series [|19| 



a s N 



71 



1 

- + 

00 



2o4iV 4 C(3) 



-4 + o(4) 



71 1 



lO 1 



00° 



Problem 16: Get eq.(81) from eq.(79), expanding it with respect to / = 7(0*). 



The first term in eq. (0) is the anomalous dimension of the GLAP equation |TTJ in leading 
order of as at 00 — > 0, which gives the solution for the structure function at x — > and 
corresponds to so called double log approximation of perturbative QCD ( DLA).The DLA 
sums the contributions of the order (aslog(l/x) log(Q 2 /<2o)) n i n the perturbative series 
of eq.(45). 



However, we would like to stress that eq. (|81|) is valid only at fixed as while the 
anomalous dimension in the GLAP equation can be calculated for running as- It means 
that we have to introduce the running as in the BFKL equation to achieve a matching 
with the GLAP equation in the region where to <C 1 and ^ < 1. 

The second remark is the fact that we can trust the series of eq. (|8T| ) only for the value 
of 00 3> ool , where 

a s N 1 4Nln2a s fatJ . 



7T 
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From eq. H7§) one can notice that moment variable N denned such that N 
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In vicinity u — > u>l we have the following expression for 7(0;) 



7M = \ + ^ u ^ 3) • (83) 

Problem 17: Get eq.(83) front eq. (79), using the property of ^(f)- 

Substituting eq. (jS3|) in eq. (|78|) we have 

^M) = / -^e^A 2 ) = [ ^e^ + ^M-i)^) = (84) 



2-Ki rv " 7 27ri 



e V * <A(o;] 



2vri 

Evaluating the above integral using saddle point approximation we obtain 

1 r 2 

us = u L + -^-—7- ■ 7-5 (85) 



^14C(3) 4y 2 



which gives the answer: 



/ n 



2 k t 



<P(y,kf) = tL.^), 2 ^"^ ,^^, . (86) 



We can trust this solution in the kinematic region where (Zn%) 2 < 2s*28C(3)y. The 
solution of eq. (|5BJ) illustrates one very important property of the BFKL equation, namely 
k\ can be not only large, but with the same probability it can also be very small. It means 
that if we started with sufficiently big value of virtuality q$ at large value of y = ln(l/x) 
due to evolution in y the value of fc 2 could be small ( k t ~ A, where A is QCD scale ). 
Therefore, the BFKL equation is basically not perturbative and the worse thing, is that 
we have not yet learned what kind of assumption about the confinement has been made 
in the BFKL equation. 

Our strategy for the further presentation is to keep A; 2 > g 2 and to study what kind of 
nonperturbative effect we can expect on including the running «g in the BFKL equation, 
as well as changing the value of lol in the series of eq. (PH). 
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10.3.4 The bootstrap property of the BFKL equation. 



We have discussed the BFKL equation for the total cross section, however this equation 
can also be proved for the amplitude at transfer momentum q 2 7^ 0, and not only for 
colorless state of two gluons in t-channel. The general form of the BFKL equation in u - 
representation looks as follows |L6| ( see eq. (|78|) : 



(co - uo G {k 2 t ) - uj G {{q-k t f))<P{uj, q ,k t ) = ^X R f^K(q,k t ,k'Mu;, q ,k' t ), (87) 

Z7T J 71 

where the kernel K(q, kt, k' t ) describes only gluon emission and 

*t (q-h) 2 q 2 



K (q, K K) = n T7T7772 + 



(h - K) 2 kf {k t -Kf(q-k' t y (ft-*O a *?' V ' 

\r is colour factor where Ai = 2A§ = N for singlet and (A^ 2 — 1) representations of colour 
SU(N) group and u) G (k 2 ) = a G (k 2 ) — 1 ( see eq. 



The bootstrap equation means that the solution of the BFKL equation for octet colour 
state of two gluons ( for colour SU(3) ), should give the reggeizied gluon with the trajectory 
a G (k 2 ) ( or oj G {k 2 ) ) given by eq. (|67D . The fact that the gluon becomes a Regge pole 
have been shown by us in the example of the next to the Born Approximation, and has 
been used to get the BFKL equation in the singlet state. It means that the solution of 
the BFKL equation in the octet state should have the form of a Regge pole : 

<P(co,q,k t ) = C0n t (89) 
us — ijj^\q l ) 

Assuming eq. Q8"9p, one arrives to the following bootstrap equation: 



u G (q 2 ) - u G (k 2 ) - u G ((q-k t ) 2 ) = ?A 8 / ^K(qMA) , (90) 



2TT J 7T 

It is easy to check that the trajectory of eq. (|67p satisfies this equation. It is interesting 



to mention that we can use eq. ( |?DD to reconstruct the form of the kernel K(q, k t , k' t ), if 
we know the expression for the trajectory ( see ref. [ 



10.4 Corect degrees of freedom at low x. 

In the previous subsections we gave the traditional derivation of the BFKL equation, which 
is rather difficult and demands a good experience in calculation of Feynman diagrams. It 
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is very important to understand better the physical meaning and the formal grounds of 
the BFKL equation. I firmly believe that during this year A. Mueller ( see also related 
papers P2[ ) has achieved the considerable progress in both understanding and formal 



derivation of the BFKL equation and its generalization. Mueller's main idea is to construct 
the small x infinite momentum partonic wavefunction of a hadron in QCD while BFKL 
calculated the amplitude for n - gluon production in so called multireggeon kinematic 
region. The wavefunction gives us much richer information on the hadron interaction, has 
very transparent physical meaning and makes the bridge between our parton approach 
to "hard" processes and new phenomena that we anticipate in the region of high density 
QCD. 



The technical trick that has been used is also very instructive, namely it turns out that 
the wavefunction looks much simpler in the mixed representation in which each parton is 
labeled by its fraction of the total hadron momentum Xj and the transverse coordinate r ti . 
The transverse coordinate is especially useful since in the low Xi region the i th gluon can 
be considered as being emitted from the system of i — 1 partons with spatial transverse 
coordinates of these "sources" being frozen during the emission of % — th gluon. Thus we 
can consider these of % — 1 partons as a system of [i — 1) qq dipoles since each gluon can 
be viewed as quark - antiquark pair if number of colours N c is big enough. So the only 
thing that one needs to write down is the emission of the i — th gluon by such a system 
of dipoles. This problem has been solved in ref.[]2T|. For example for emission of gluon 



(#2, r tc = r 2) from one dipole which is the quark (x q — 1 — x%, r^q — To = 0) and antiquark 
(xq = Xx, r t q — ri) is equal to 

^(x 1 ,x 2 ;r 1 ,r 2 ) = -^ (0) (*i; rr){^ - ^} - e\ . (91) 



7T 
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2. 

Mueller made one very important step in our understanding of our parton system, 
namely he found what sum rule plays the role of the momentum sum rules in the GLAP 
approach for low Xj partons. This sum rules is the normalization of the partonic wave- 
function: 

JW 1 ^\{d 2 r u \^{x 1 ,...x n -T tl ,...T tn )\ 2 = 1. (92) 

Using this equation one can easily take into account so called virtual corrections, which 
in this case are mostly known as gluon reggeization or non -Sudakov form factor. The 
importance of this step can be compared only with the transition from Gribov - Lipatov 



43 



form of usual evolution equation with Sudakov form factor in the kernel to well known 
Lipatov - Altarelli - Parisi elegant form based on direct use of the momentum sum rules 
in QCD. 

3. 

The physical application of this new approach has not been considered but Mueller 
noted at his Durham talk |23| that his approach will be able to resolve the old problem 



with the BFLK equation. Indeed, the physical meaning of the growth of the structure 
function at Xb — > is the increase of the number of "wee" partons ( N ) that can interact 



with the target ( N oc x B u ° )( see, for example, review ||24|| ). However the multiplicity of 
gluons calculated as the ratio / J^ 7 . J &t turns out to be small ( of the order of a s ln^ 
). It means that partons are in a very coherent state in a typical inelastic event. However 
the behaviour of the parton cascade at large multiplicity should be quite different from 
Poisson distribution since at large multiplicity all N parton can be freed in the interaction. 



The main progress was related to the fact that they discovered the correct degrees 
of freedom for the region of small x in QCD. They showed that if we discuss the deeply 
inelastic processes not in terms of quark and gluons as we did before but introducing 
new degrees of freedom, namely the color dipole of the definite size (of), we still can use 
the simple probabilistic interpretation in the region of low x. It means that the physical 
meaning of the deep inelastic structure function is not the number of quark or gluons but 
the number of the color dipoles with the sizes larger than 1/Q (d > 4). It is worthwhile 
mentioning that at x ~ 1 we still can use our old interpretation as probability to find 
parton but as well as a new one. 



Problem 18: Read Mueller's paper[21] and try to understand how simple looks the 



BFKL equation in new degrees of freedom. 



The main progress in the Mueller approach was related to the fact that he discovered 
the correct degrees of freedom for the region of small x in QCD. He showed that if we 
discuss the deeply inelastic processes not in terms of quark and gluons as we did before 
but introducing new degrees of freedom, namely the colour dipole of the definite size (d), 
we still can use the simple probabilistic interpretation in the region of low x. It means 
that the physical meaning of the deep inelastic structure function is not the number of 
quark or gluons but the number of the colour dipoles with the sizes larger than 1/Q 
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(d > 4). It is worthwhile mentioning that at x ~ 1 we still can use our old interpretation 
as probability to find parton but as well as a new one. 

This idea shed a light on all technicalities of the BFKL equation which was the mani- 
festation of the artistic skill in the calculation of Feynman diagrams and reduced the high 
art to normal physics understandable for a normal post-doc. I anticipate a big progress 
in the development of all kind problems related to high energy QCD based on Mueller 
approach. 



11 New physical phenomena at small xb- 

Let us outline the new phenomena which we anticipate to occur in the region of small 
xb in perturbative QCD for the case of deeply inelastic scattering. Three of them are 
particularly important and determine the physical picture of the parton evolution, or 
cascade as we shall call it, in the region of small xb, namely: 



1. The increase of the parton density [16| at xb — > 0. 



2. The growth of the mean transverse momentum of a parton inside the parton cascade 



at low 16 25 . 



3. The saturation of the parton density 

Before discussing each of these phenomena in turn, let us try to understand these new 
phenomena at small xb by recalling some facts regarding perturbation theory and deeply 
inelastic scattering. According to the factorization theorem, the deep-inelastic structure 
function is a product of a parton density and a short distance coefficient function, which 
is calculable in perturbation theory. Both factors depend on the factorization scale, which 
we choose equal to Q 2 . Concentrating on the gluon density, in the kinematical region we 
are considering the former can be represented as 

xG(x,Q 2 ) = ^ n C n (Ql)(a s L) n + 0(a s (a s L) n ) , (93) 

where L is the large logarithm in our problem. The coefficients C n (Qg) contain nonper- 
turbative information and depend on the initial scale Qq of the cascade, see Fig. 28, 

To see what kind of large logarithm could occur, one can examine the probability of 
emission Pj of the i-th parton in the cascade. Near x = it can be written as 



Pi 



N c a s dxi dqf t 

7T Xi qf t 
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1. In the region of large virtualities of the photon in deep inelastic scattering (but not 
at very small values of xb) Lis equal to InQ 2 . 

2. At Q 2 fixed (Q 2 ~ Ql 3> rn 2 roton ) but at small value of xb we have a different 
large logarithm in eq.(2), namely L = ln^-. 

3. When both InQ 2 and In are large (Q 2 ^> Ql xb — *• 0) the large logarithm has a 
more complicated form, namely L = InQ 2 In In this case one can apply the so-called 
double logarithmic approximation (DLA) of pQCD. The objective now is to resum these 
large logarithms in eq.(9), so that we can estimate the behavior of the parton densities 
in the corresponding kinematical regimes. This resummation is performed via evolution 
equations, which differ depending on the type of large logarithm under consideration. 

11.1 Increase of the gluon (quark) density at xb — > 0. 

We now investigate the new phenomena at small xb, but begin by discussing the evolution 
equations that sum particular classes of large logarithms in the parton densities. 

Gribov - Lipatov - Altarelli - Parisi ( GLAP) evolution equation. 

The GLAP equation sums those contributions in the parton cascade which compensate 
the smallness of the QCD coupling constant a s by the large logarithm InQ 2 . It reads 



where P(z) is the Altarelli-Parisi evolution kernel, calculable in perturbation theory, and 
<p(x,Q 2 ) is a generic parton density (matrix). 

Clearly, from the expression for P h we need strong ordering (Q 2 » ... >■ q 2 t >> 3> 

-jp = Qq) in transverse momenta of emitted partons since only under such a condition does 

one get a In Q 2 contribution for each integration over q it in the parton cascade in Fig. 28. 
This strong ordering means that the GLAP evolution equation allows us to calculate the 
probability to find a parton of transverse size r t ~ -k inside the initial parton (quark or 
gluon) in the hadron at fixed xb (see Fig. 29). In the case where both Q 2 and 1/xb are 
very large (the DLA approximation), the AP kernel simplifies to P(z) = N c a s /irz, and 
the solution to the GLAP equation for the gluon density reads then 



d 



<P(x,Q 2 )=a s (Q 2 ) f ld ±P{ z )^,Q 2 ) 

Jx Z Z 



(94) 



dlnQ 2 



xG(x,Q 2 ) ~ xG(lo ,QI) exp( 
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where ui = In ^/ In \ ) 1 ^ 2 , and Q is the initial scale of the parton cascade. 



Problem 19: Find the exact solution of the GLAP equation at low x ( see eq, (94)), 
using P(z)Uo - ^f- 



Note that the solution of the GLAP equation decreases as Q 2 increases. This agrees 
with the intuitive notion that, because of asymptotic freedom, the hadron should be 
almost empty at small distances. 

Balitski- Fadin -Kuraev - Lipatov (BFKL) evolution. 

The situation changes crucially if we consider the behavior of the parton density <ft 
in another extreme regime: at fixed Q 2 but small xb- Here we have the large logarithm 
L — In — . The picture of the parton cascade, Fig.30, shows that the number of partons 
increases drastically in the region of small xb, since each parton in the basic branching 
process in Fig. 28 is allowed, due to the abundance of available energy, to decay in its own 
chain of daughter partons. Let us estimate the total multiplicity of gluons Nq associated 
with this complicated process. 

We are able to estimate the number of cells in this chain diagram, which corresponds 
to the typical number of parton (gluon) emissions in one snfrprocess in our parton cascade 
(jTg). The characteristic value of Ay iji+1 , the rapidity difference between two adjacent 
rungs of the ladder, is equal to Ay iji+1 oc l/a s , from the expression for gluon emission 
Pi. Thus 

_ ln^- l 

n G = XB oc a s ■ In — . (95) 
Ayi,i+i %B 

The total number of Nq partons that could interact with the target can then be calculated 
as follows (see Fig.30): 

N G oc e 7 ^ = e caslnl/xB = (— ) CQs , (96) 

x B 



where the constant c should be calculated using the exact BFKL equation [0. In first 
approximation ca s = 0.5. 
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11.2 A new scale for the deeply inelastic process at xb — > 0. 

Thus in the region of small xb the density of partons increases. A careful study of the 
behavior of the parton cascade shows that also the mean transverse momentum of the 
parton increases. The reason is as follows. In the case where we neglect the running 
coupling constant of QCD our theory is dimensionless so each emission leads to a value of 
the transverse momenta of the daughter gluons which is of the same order as the transverse 
momentum of the parent gluon. This could be seen just from the explicit expression of 
Pj. We introduce the ratio \\nqf t /q 2 +lt \ which characterizes the emission. It is roughly 
constant. 

After no emissions 




< \HTn.t/Ql)\ > oc Vn G ~ Wot, In— , (97) 



since the parton cascade corresponds to a random walk in the variable ln(g„ t /Qo) BUI [26 
(see Fig. 30). Although we will not explicitly present the BFKL equation here, we will give 
the solution to this equation with the initial condition xbG(xb = %o) = 8(\n(Q 2 /Qq)). 
This is instructive because it shows explicitly all the properties mentioned in the above: 

G(y -yo,r- r ) = \j ^ * / 1 = * ex P{Mv ~ Vo) - g x~ 7 ^ — 




where 

a s (Q 2 o)N c ^ 14C(3) x Q 2 

u = 4 In 2, A = , y-y = \n — , r - r = In — ^ 

7r 4 In 2 x B Qq 



Problem 20: Calculate the mean value of \n(Q 2 /Ql) using the solution of the BFKL 
equation. 



11.3 Saturation of the gluon density. 

The increase of the parton density leads to a new problem in deeply inelastic scattering, 
namely the violation of s-channel unitarity. This is the requirement that the total cross 
section for virtual photon absorption be smaller than the size of a hadron. 

a( 7 *iV) < txR\ (98) 
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At small xb the gluon density dominates and we can write 



ah*N) = const^- xG(x, Q 2 ) (99) 
Q 

where const a em /Q 2 approximates the cross section cr(7*g). We have shown previously 
that the value of the gluon density increases very rapidly as x — > 0. Using the DLA result 
(see eq.(ll)), we can rewrite the unitarity constraint in the following way, dropping the 
const: 

Q 

Q 2 y x B ' 

Here we replaced a em by a s since the probe can also be a virtual gluon, not just a virtual 
photon. 



• < nR%. (100) 



From this expression alone one can conclude that unitarity will be violated PSl at 



x < x cr where log - = - • — H^j) (101) 

Therefore unitarity is violated even for (very) large values of Q 2 when x < x cr . Clearly 
the miraculous confinement force cannot prevent this from happening. Thus we have to 
look for the origin and solution of this problem within pQCD. 

Let us try to understand what happens in the region of small Xb by examining the 
parton distribution in the transverse plane (see Fig.29) Q. Our probe (photon) feels those 
partons whose size is of the order of 4. To begin, at x ~ 1 we have only several partons 
that are distributed in the hadronic disc. If we choose Q 2 such that 

« Rl (102) 

then the distance between partons in the transverse plane is much larger than their size 
and we can neglect the interaction between partons. The only process which is essential 
here is the emission of partons that is taken into account in the usual evolution equation. 
As x decreases the number of partons increases and at some value of x = x CT partons 
start to densely populate the whole hadron disc. For x < x cr the partons overlap spatially 
and begin to interact throughout the disc. For such small x-values the processes of 
recombination and annihilation of partons should be as essential as their emission. Both 
these processes are however not incorporated into either the GLAP or BFKL evolution 
equations. 

^Recall that a high energy hadron in the parton model can be represented as a Lorentz contracted 
disc. 
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What happens in the kinematical region x < x cr is anybody's guess but there is 
enough experience with some models [27] to suggest the so-called saturation of the 
parton density []25|| . This means that the parton density 0, which we discussed before, 
is constant in this domain. 



12 Shadowing correction and nonlinear evolution equa- 
tion. 

12.1 Qualitative Derivation of Nonlinear Gribov-Levin-Ryskin 
Evolution Equation 

In section 11.3 we remarked that at small x (high energies), the density of gluons becomes 
so large that the unitarity constraint is violated, even at large values of Q 2 . We argued 
that the physical processes of interaction and recombination of partons, which are usually 
omitted in perturbative calculations, become important in the parton cascade at a large 
value of the parton density. To take interaction and recombination into account we must 
identify a new small parameter that lets us estimate the accuracy of our calculations. 
This small parameter is [^ 



O.s 



W = ^p(x,Q 2 ) . (103) 
The first factor is the cross section for absorption of a gluon by a parton in the hadron; 



the transverse density p is defined in eq. ( 54 ) . Effectively, W is the probability of a parton 



recombination during the cascade. We can rewrite the unitarity constraint eq. (|100|) in the 
form 

W < 1 . (104) 

Thus, W is the natural small parameter in our problem. Amplitudes that take gluon 
recombination into account can be expressed as a perturbation series in this parameter 
2l|. We can resum this series, and the result |25j can be understood easily by considering 



the structure of the QCD cascade in a fast hadron. Two processes occur inside the cascade 

emission (1 — > 2); probability cx « s| o ; (105) 



annihilation (2 — > 1); probability oc a 2 s r 2 p 2 oc a 2 — p 2 

Q 

where r is the size of a parton produced in the annihilation process. For deep inelastic 
scattering r 2 oc l/Q 2 - 
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At x ~ 1 only the production of new partons (emission) is essential because /> C 1, 
but at x — > the value of p becomes so large that the annihilation of partons becomes 
important. This simple parton picture allows us to write an equation for the density of 
partons that properly accounts for all these processes. The number of partons in a phase 
space cell (Ay = Aln(l/a;), AlnQ 2 ) increases through emission and decreases through 
annihilation. Thus the balance equation reads 

d 2 p a s N c a 2 7 



p - -krp\ (106) 



din ± din Q 2 7T H Q 2 

or in terms of the gluon structure function xG(x, Q 2 ) 

d 2 xG{x,Q 2 ) a s N c 2 a 2 7 , nt ^2^2 / 1fV7 N 

din 1 ding 2 = —vG^Q ) " ^OHP^'Q )) • ^ 107 ) 

This is the so-called Gribov-Levin-Ryskin (GLR) equation p5| . The parameter 7 can in 



fact be calculated order by order in W perturbation theory. Mueller & Qiu |28| found it 
to be 

7 = — foriV c =3. 
' 16 



12.2 The Scale of the Shadowing Corrections 



The second term in Equation |107| describes shadowing corrections (SC). Its size crucially 



depends on the value of R 2 . The physical meaning of R 2 is clear: R is the correlation 
length between two gluons in a typical hadronic situation (at x ~ 1). In our derivation 
of the GLR equation we assumed that there are no correlations between gluons except 
those imposed by their confinement in a disc of radius R. 

If R = Rh the value of the SC is negligibly small. 



If R -C Rh, the SC could be large (see Refs. p9[] , |30| and [[H]]). Recently Braun et 



al |32j performed the first theoretical estimates of the value of R within the framework of 
QCD sum rules. They found 

R = 0.3 -0.35 fm ~ - R h . 

3 

an encouraging result for experimental study of SC at present energies. For the gluon 
structure function the solution of the GLR equation is shown in Fig. 32 and the value 
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of xG(x, Q 2 ) from which saturation starts depends on the value of R 2 . For R = R p 
saturation starts from 

xG{x,Q 2 ) > 150 at Q 2 = 10GeV~ 2 

while for R= ^R p 

xG(x,Q 2 ) > 15 at Q 2 = lOGeV- 2 . 

It should be stressed that the GLR equation predicts the value of the maximal packing 
factor ( PF) which is 

PF =< \r 2 constituent \ > -p (108) 

PF does not depend on the value of the radius R and for parton with < \r 2 onstituent \ > = 
^2 the GLR equation gives: 

(PF) max = 0.21 for N c = 3 



12.3 New Scale of Transverse Momentum (the Critical Line) 

12.3.1 The main property of a solution of the GLR equation 

First, let us rewrite the GLR equation (see eq.( |107| )) in terms of new variables: F = 
xG(x,Q 2 ); y = SN c /b ■ ln(l/x) and f = lnln(Q 2 /A 2 ), as 

^£ = - 7'exp[-e - e«]F} , (109) 

dydt, 2 

where 7' = 2 717 jbR 2 A 2 . The terms ( — £) in the exponential correspond l/Q 2 (at s ). 
Note that the term proportional to 7' corresponds to the recombination probability W, 
defined below Equation (26). By inspecting eq.( |109| ) we can guess the main property of the 
solution. At the start of the evolution, when the value of the structure function is small, 
we can neglect the nonlinear term. The derivative with respect to y is here positive, so 
the structure function increases monotonically with y, or ln(l/a;). This growth continues 
up to the value 

Fmax = — ■ exp[i + e ? ], 

7 

where the derivative vanishes. If F > F max the derivative from the GLR equation becomes 
negative and F max is thus the limiting value of the gluon density and the solution of the 
GLR equation as x — > 0. 
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Therefore, at each value of Q 2 there is a critical value of x = x cr when F = F max or, 
conversely, at each value of x there is a critical value of Q 2 = q^ (x) when F = F max . This 
value Q 2 = q^x) introduces a new transverse momentum scale in the parton cascade. 
The behavior of the gluon density changes crucially at Q 2 = 0{cfo{x)). Indeed, 

F = F GLAP (x, InQ 2 ) for Q 2 > q 2 (x), 

where F GLAP means that F is a solution of the linear GLAP equation and a smooth 
function of InQ 2 . Alternatively 

Q 2 
Qo\ x ) 

We postpone the discussion of the physical meaning of this new scale and instead deter- 
mine first the x dependence of q^x). 



12.3.2 Trajectories 

A semiclassical analysis offers new insights into solutions of the GLR equation. It is easiest 
to find the solution to the linear equation by passing to a moment (cu) representation and 
introducing the anomalous dimension 7. In moment representation the solution is as 
follows: 

F(x,Q 2 ) = /^Mo^QV^ 7 ^ , (HO) 

where M (u,Q 2 ) is the initial distribution at Q 2 = Qq and the anomalous dimension 
7(o>) = \/2uj. At large y (small x) we can evaluate the integral in the saddle point 
approximation. The position of the saddle point is determined by: 

= d\nM (u) _ rfTM ^ (m) 
duo duo 

This equation describes a family of curves y = y(£) that can be considered as semiclas- 
sical trajectories of parton evolution because u is a constant of a motion in the linear 
evolution equation. The semiclassical values of F along the trajectory give an approxi- 
mate estimate of the structure function at fixed Q 2 and x. Although this method yields 
no new understanding for the linear equation, it can be generalized to the nonlinear case. 
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12.3.3 Semiclassical Solution to the GLR Equation 



In semiclassical approximation we parametrize the solution to the GLR equation by writ- 
ing 

F = e , assuming S y ^> S y ^ , 

where S y = dS/dy , S% = dS/d^ and S y ^ = d 2 S/dyd£. The solution to the linear 
equation is S = y/2y£,, for which one can check that the above inequality holds to good 
accuracy 



Thus in the semiclassical case the GLR equation can be reduced to the form 

1 
2 

Eq.( |112|) can be solved using the method of characteristics, or trajectories (|25j |33 



S y Sz = -(1 - V exp{S - e ? - £}) • (H2) 



p4|). By introducing an intrinsic time coordinate t along the trajectory and denoting the 
derivatives with respect to this time variable by y etc we can rewrite the GLR equation 
as the following system of ordinary equations: 

V = S S 

Z = Sy 

S — — 2 Sy 

S y = — 7' Sy exp{ S — — £ } 
St = -y [5 C - 1 - e ? ] exp{S - e ? - . (113) 

In the linear case (7' = 0) this system of equations has a simple solution: 

S v = oj = Const ; Se = — = Const , 



which coincides with eq. ( |111| ). We introduce a; here such that it corresponds to u in eq. 
(111). From eq.(113 ) we find 

y - y(t = 0) = ±-t , £ - £{t = 0) = cot , S = h 

Thus we obtain 

S=^2(y-y(t = 0M-at = 0)) , y - y(t = 0) = -L(£ - £(t = 0)) 

The second of these equations corresponds to Equation 33 with 7 = 1/cu, the correct 
anomalous dimension in the GLAP equation for x — > 0. Along each trajectory of the 
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linear equation, the rescattering probability W = 7' exp{S — — £ } increases at the 
beginning when S ~ 1, reaches a maximum value when S% = (1 + ) and then decreases 
at larger £. However, the trajectories of the GLR equation become quite different from 
the linear ones. Indeed, in the region where W becomes of order one, the derivatives S y 
and S% are not small and trajectories are not straight lines. A number of papers ( ||25|| EJ 
[Q) have shown that critical line separates the region of almost-linear trajectories (small 
values of W < 0(a s )) and strongly nonlinear solutions (see Fig. 33). 



It is easy to check explicitly that for large £, 

S = e $ - ln 7 '; S y = e"*; « ^ (114) 



is the solution of eq. ( |113j ) and of the GLR equation as well. To the right of the critical 



line (x > x cr ) the equation y cr = ^- may be recast as 

= m hW/A2) - (115) 

The nonlinear corrections are so small that we can neglect them and use the trajectory 
of the standard GLAP equation. However, when x becomes smaller than x cr , W — > 1 
and we have to use the nonlinear trajectory. Remarkably, this trajectory cannot cross the 
critical trajectory (see Fig. 33). Using this property, one can suggest the following way to 
solve the GLR equation. To the right of the critical line it suffices to find the solution of 
the linear (GLAP) equation with the new boundary condition F = e s with S = e^ — lnj' 
on the critical line (see eq.( |114] )). To the left of the critical line we have a separate system 
of trajectories. The solutions to the left do not depend on the solutions to the right of 
the critical line. Here we enter the very interesting kinematic region of hdQCD. Even the 
GLR equation is not the right tool to solve problems in this region because in addition 
to the recombination of two partons, simultaneous interactions of three, four and more 
partons are also important in this region. 



12.4 The Physical Meaning of the New Scale for the Typical 
Transverse Momentum in the Parton Cascade 

The equation for the critical line (eq.( |115| )) introduces a new scale for the value of the 
typical transverse momentum in the parton cascade, namely 

q \ = g 2 (x)|_ - A^eVW^. (n 6 ) 
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This new scale plays the role of an infrared cutoff in inclusive production. This new 
infrared cutoff has a dynamical origin: parton-parton recombination in the parton cascade. 
It is not related to confinement. 



The quantity qo(x) is also the Landau- Pomeranchuk momentum [[36] in the parton 
medium. Indeed, in the saturation region q t < qo{x) the mean free path A does not 
depend on x. In this region A = 1/aT where a is the cross-section of parton-parton 
interaction (cr oc a s l/q 2 ) and T = xG(x,q t ) is the number of partons with transverse 
momentum q t at fixed x. If q t < qo{%) T oc q 2 <po and oT = Const(q t ). Thus A is equal 
to the formation length of a parton with x = xq, which can be derived from the relation 
qo(xo) = Qt- If qt > Qo(x) then aT oc xG(x,q t )/q 2 <C 1/mx so A is big. Therefore, qo(x) 
is precisely the LP momentum, and the factor C that interpolates between coherent and 
incoherent emission of gluon with transverse momentum q t and fraction of energy Xj can 
be rewritten as follows: 

C = ^ p^r . (117) 

qt 



1 + 



12.5 Correlations 

We assumed that the only correlations between gluons in the cascade arise from the fact 
that the gluons are distributed in a disc of radius R. The probability to find two partons 
with the same value of ln(l/a;) and InQ 2 is 

P 2 oc p 2 . (118) 

However this assumption is only correct if the number of colors is large (N c — > 00) [|35|j . 
If N c is not large, Bartels and Levin et al claim |35j claim that 



P2 ( MiV^ 1 Q 2 1 

— -> exp (4/ In - In -2 — i ) . (119) 



Bartels & Ryskin |37| found that this correlation reduces the correlation radius R by a 



factor of 1.3 - 1.4 in the HERA kinematic region. 

It is also very instructive to write down the evolution equation that takes into account 
the correlation of eq.( |119| ). This is a system of two equations, instead of the single 
equation of eq. ( |107| ), namely 

92 P asNc ~ a ^P 2 , (120) 



dln^dlnQ 2 vr F Q 2 

d 2 Po 4a N 1 2a 2 ^ 

" c (l+ * t j P a ~ -^PiP- (121) 



d\n±d\nQ 2 7r v (N 2 -1) J z Q" 2 
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Here we assume that the probability to find three partons in the same cell of the parton 
phase space is P 3 = P 2 • p. Without this assumption the last term in Equation ( |121| ) reads 
as 2a 2 j/Q 2 ■ P 3 . The general equation for P n is (see Ref.pTf for details) 



d 2 P n na?7 

This equation indicates the need for the solution of two theoretical problems. 

(i) We should find a general equation for 7 n , which has the formal meaning of the 
anomalous dimension of a high twist (2n) operator (Refs. [pT| ,|38||, [j35H ). 

(ii) We need a closed form of the equation which incorporates all P n for the deep- 
inelastic structure function, and its solution. 



The first problem was solved by A.G. Shuvaev and us in Ref.pH]. The complicated 



problem of gluon-gluon interaction was reduced to that of interaction of colorless "gluon 



ladders" (Pomerons). It was shown in Ref. ]35| that this idea works for the case of the 



anomalous dimension of the twist four operator, in other words, for the case of P 2 . The 
fact that we can consider the rescattering of n Pomerons to find 7„ (or P n ) signifies that 
we are dealing with a quantum problem: the calculation of the ground state energy for an 
n-particle bosonic system where the interaction is attractive and given by a four-Pomeron 
contact term, with a coupling proportional to A = 4a s / (N 2 — 1). This observation simpli- 
fies the problem and enabled us to reduce it to solving the nonlinear Schrodinger equation 
for n Pomerons in the t-channel. It should be stressed that this effective theory is a two- 
dimensional one, the two dimensions being spanned by ln(l/x) and In Q 2 : the Schrodinger 
equation can be written for n particles on a line. It is well-known that this problem can 
be solved exactly (e.g. by Bethe Ansatz) and the answer was found to be 

N ca s n 2 1 1 2 

con s (iVp — \y 

Although Pomerons are bosonic, maximally two Pomerons occupy a single-particle level 
in this groundstate configuration, thus behaving like fermions. This is due to the two- 
dimensional nature of this problem. The second problem has also been (partly) solved by 
us (unpublished but see Ref.[|3ll for results). We introduce the generating function 



P(x,Q 2 ,r ] ) = J2Pne nv (124) 



n=l 



with Pi = p. Using eq.(fL22|) we obtain 



«91n±«91nQ 2 vr ^ + 3(^2.1)2^™ P 



-o^e-MQV^-P), (125) 
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where P v = dP/dr], P m = d 2 P/drf etc. The deep-inelastic structure function is the 
solution of this equation at 77 = — ln(Q 2 ), viz. 

x G(x, Q 2 ) = ttR 2 P{x, Q 2 ,v = ~ HQ 2 )- (126) 

Eq. ( |125|) is a new evolution equation which allows one to take into account both initial and 
dynamical correlations between gluons in the parton cascade. The key problem in solving 
it is nonperturbative in nature, namely to determine an initial condition at x = x°. The 
dependence of P(x°,Q 2 ,r]) on 7] describes the correlations between gluons in the hadron 
near x ~ 1. The simplest case is to neglect this initial correlation and use the eikonal 
approach for P(x°,Q 2 ,r]), namely 



P(x°,Q 2 , V ) = J2e n ^—^P 1 (x°,Q 2 ) n = 1 -exp(-e"P 1 (x ,Q 2 )) (127) 



The solution of this equation has not yet been found (see Ref.[|3T] for more details) 



13 A look at the first small x HERA data. 
13.1 25 nb of Experimental Data 

Let us list conclusions from the first experimental results at small x from the HI and 
ZEUS collaborations at HERA (||): 

1. The increase of the value of the deep inelastic structure function F 2 at x — > 1CT 4 
shows that the gluon density is large: 

xG(x,Q 2 ) -> 40 - 50 at x B = 10" 4 and Q 2 = 20GeV 2 . 



The value for xG(x,Q 2 ) ~ 50 was determined from the MRSD-' [40] parameterization 
that roughly describes the new data. 

2. The ZEUS collaboration measured the diffraction dissociation (DD) f\ cross section 
and found 

= ^2nb = g x 10 -2 atx = 1Q -4 and Q 2 > 1QGeV 2 
a-y* p 80 nb 



"DD is the process in which a number of hadrons are produced in association with a recoil proton, 
whose energy fraction is close to one. 
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where <7 7 . p is the total cross-section for virtual-photon absorption. We will say more about 
the diffractive component shortly. 

3. The structure function F 2 can be described by the GLAP evolution equation as- 
suming that at the initial virtuality Qq ~ 5GeV 2 the gluon structure function increases 
rapidly at x — > as xG(x, Q%) oc x~ x and A > 0.3. Note that at sufficiently large Ql A 
approaches uq, which occurs in the solution of the BFKL equation (see section 11.2). 

4. Both collaborations measured the total photoproduction cross section at for very 
small virtualities of the incoming photon and also for real photons. This cross section 
exhibits a common property of soft processes, namely a t — > s A where A ~ 0.08. 



13.2 The Value of the Shadowing Corrections from HERA Data 

Let us first try to understand the large number of gluons (between 40 and 50 at x = 10~ 4 )). 
An interesting way to do this is to compare the proton at x — 10~ 4 and Q 2 = 15GeV 2 
with an iron nucleus, which has a similar number of nucleon constituents. If we calculate 
the packing fraction by assuming the size of the constituents to be Rh = 0.86 fm for 
the nucleons and ra = 1/Q for the gluons, we find that the two packing factors are very 
similar if one assumes the size of the gluon disc to be -R/i/3 (see section 4.2). Thus, in this 
sense, we see that a proton x = 10~ 4 and Q 2 = 15GeV 2 is very similar to an iron nucleus. 

We now estimate the size of shadowing correction AF 2 for the deep inelastic structure 
function. Let us write it as 

F 2 (x,Q 2 ) = F^ LAP {x,Q 2 ) - AF 2 (x,Q 2 ) , (128) 

where F 2 GLAP is the solution of the usual (GLAP) evolution equation. A straightforward 
estimate gives 

— ^ oc a s < r 2 G > p -> 0.1 for Q 2 = 15 GeV 2 and x = 10~ 4 . (129) 
F 2 

However, the ZEUS data on diffraction dissociation let us better estimate the value of 



SC. Indeed, Levin & Wuesthoff [41] showed that 



|AF 2 | = F 2 DD , (130) 

where F 2 DD relates to aER in the same way that F 2 relates to a r * p . In other words 
|A<7 7 * p /cr 7 *p| = c££ , /(T 7 *p. From the ZEUS data we can then conclude that | AF 2 \ / F 2 > 
6 ■ 10~ 2 . Note that in this case we can begin not only to discuss whether there are SC, 
but even whether we could predict this SC value in our theory. 

The details of the evolution equation for F 2 DD were discussed in Ref. jJTJ, in which 



Levin & Wuesthoff suggested measuring the sum F 2 + F 2 DD in which all contributions of 
SC cancel. One can use the GLAP equation for this sum, even at small x. 
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13.3 Saturation of the Parton Density or Different Physics for 
"Soft" and "Hard" Processes? 



The HERA data allow us to raise the question that is the title of this subsection. Indeed, 
we can distinguish two regions in Q 2 that have quite different energy behavior of the 
photoproduction cross section. The experimental facts that pertain to these regions are 
as follows: 

1. Q 2 <C lGeV 2 . (a) The total cross section is essentially constant here (or in- 
creases slightly with energy); (b) <j DD /a t ~ 10%. Thus in this kinematic region the 
photoproduction process seems to be a typical soft process like a hadron-hadron collision. 

2. Q 2 > 5GeV 2 . (a) The total cross section increases rapidly with energy (x) 

A ; (6) cr^,p/cr 7 .p ~ 10%. This is the typical deeply inelastic scattering process. 
Both experimental facts arise from the contribution of a so-called hard Pomeron, which 



is a solution to the BFKL equation [TIJ and leads to o" 7 . p oc x u ° with lu > 0.4 while 
the smallness of the ratio cr5^/er 7 * p has a natural explanation in the small values of SC. 

What happens for intermediate value of photon virtualities lGeV 2 < Q 2 < 5GeV 2 7 
There are two different scenarios for this kinematic region: 

1. Landshoff picture At Q 2 < lGeV 2 all experimental data for photoproduction 
as well as for other soft processes can be described by the exchange of the soft Pomeron, 
which is a standard Regge pole, with intercept ap = 1 + e (e ~ 0.08 <C u ) (see Ref.H). 
The small value of the ratio o DD jo t can then be interpreted as an indication that the 
SC are small and can be treated perturbatively. In particular Donnachie and Landshoff 
considered only two Pomeron exchanges. In this scenario the region of deeply inelastic 
scattering has different underlying physics described instead by a hard Pomeron. The dis- 
tinction between the hard and soft regions seems quite arbitrary at present. However, this 
approach is very simple and provides an elegant description of all available experimental 
data on soft processes. 

2. Saturation of the parton density The second scenario is intimately related to 
the hypothesis that the parton density saturates. In this scenario the hard Pomeron is 
responsible for the behavior of the total cross-section in both kinematic regions, but the 
small value of the ratio a DD /a t is interpreted differently for large and small virtualities: 
For large virtualities it supports a small value of SC, but at small ones we interpret it as 
an indication that the SC become very large and lead to a black-disc constituent quark. 
The greatest advantages of this scenario are the unified description of both the diffractive 
and inclusive processes, and that its solid theoretical background, namely properties of 
the hard Pomeron and the shadowing correction in perturbative QCD. It is worthwhile to 
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mention that strong SC correction gives rise sumalteneously the saturation of the parton 
density, smooth behavior of the total cross section at high energy, a natural explanation 
of the transition from steep energy behavior of the "hard" Pomeron to smooth energy 
dependance of the "soft" total cross section and small value of the diffraction dissociation 
cross section in both kinematic regions (see ref. [Q for more information ). However the 
data have not yet been described within this hypothesis. 

These two scenarios result in different behavior in the region of intermediate virtualities 
of photon. In the second scenario we expect some transition region with smooth behavior 
of a t vs Q 2 . A first attempt to extract this behavior from available experimental data 



shows that such a transition region does not contradict such scenario P2"| , but it is still 
too early to draw a definite conclusion. 



14 How to measure the high density event. 

Let me list here the main ideas how to measure the new physics that we anticipate at 
high density system of partons: 

1. The probability of double parton interaction should be large (of the order of the 
maximum value of the packing factor ) about 20%. 

The double parton interaction can be seen not only as cross section for production 
of two pair of hard jets with the same value of rapiditybut also as a cross section of 
the inclusive production of hadrons in the window of rapidity y + Ay, y — Ay where y 
is the rapidity of a hard jet with transverse momentum p t and Ay = hi^, where po 
is the transverse momentum of produced hadron. It could be also seen as a long range 
correlation in rapidity between produced hard jet and produced hadron which is not 
specially hard. 

2. In the high density event we should see the Landau - Pomeranchuk suppression 
of the emission of gluons with transverse momentum smaller that the typical momentum 
<1o(xb) which can be found from the equation: 

xG(q 2 (x) } x) 

qUx)*W ~ {PF)maX ■ (131) 

Such a suppression can be seen as the deviation from the factorization theorem for jets 
with p t < q (x). 

3. . Decorrelation effect for jets with transverse momentum p t of the order of go- The 
value of transverse momentum for such a jet is compensated not by one jet in the opposite 
direction but by a number of jets with average transverse momentum about go- 
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4. Polarization of produced hadrons allows us to measure the typical transverse mo- 
mentum in the process since in the region of pQCD polarization should be equal to zero. 

5. It is seen directly from eq.(9) that the saturation reaches in the system with small 
size at larger transverse momentum (smaller value of the gluon structure function). So 
this is why we have to create experimentally such compact system. We have three ideas 
how to confine the gluons in the disc of the small size: 

A. to find a carrier of partons with small size. Even hadron could be such a carrier if 
the hypothesis of constituent quarks with small radius will be confirmed experimentally. 
However better to use the virtual photon or Pomeron. The last is not well theoretically 
defined object and what is Pomeron is one of the questions that the future experiment 
should answer. However even available experimental information confirm the idea that 
Pomeron's size is much smaller that hadron one and of the order of Rp ~ \J~olp ~ 
Q.bGeV~ l ~ O.lFm. Thus hard diffraction with Pomeron can give a good possibility to 
localize the parton system in small disc and to see high density phenomena in the most 
clear way. 

B. to find the experiment (microscope) that can resolve the small part of the hadron 
and investigate it in detail. This idea is realized in so called Mueller - Navalet process 



[g§] or "hot spot" hunting |45| (see Fig.34) 



This process allows us to measure the small cx — part of the hadron and using the 
two jets with sufficiently large transverse momentum as a trigger e can study the system 
with large parton density in many details. 



C. Bjorken |46[ pointed out that the large rapidity gap ( LRG) processes can give 
us new way to look inside the high density parton system. Indeed, due to intimate 
relation between inelastic processes and elastic one coming from new reggeon-like approach 
the process with the LRG such as two high transverse momentum jet production with 
rapidities yi and 1/2 but without any hadron with rapidities y\ > yh > 2/2 can be described 
as the exchange of "hard" Pomeron (see Fig. 35). The properties of the " hard" Pomeron 
exchange is well known theoretically (see ref. [TIB] ) and can be checked experimentally. 

6. One of the way to measure the high parton density is to select the event with large 
multiplicity of produced hadrons. More discussion on this subject you can find in the 
report of Fermilab Working Group on QCD [[17} 
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TOMORROW: 



EFFECTIVE THEORY at HIGH 

ENERGY 
from NONPERTURBATIVE QCD 

15 Summary. 

Let us summarize where we stand in the theoretical development at high energy physics 
and what of the fundamental problems listed in section 9.3 have been solved. 

As we have discussed in section 9.5 we struggle with two major problems in pQCD: 
(a) Our perturbative series have a natural small parameter as, but its smallness is often 
compensated by large logs, which are different in different processes. In the case of the 
processes at high energies ( small x) the real parameter of the perturbation series is 
CKglii(l/a;) which is not small, (b) The perturbation series is asymptotic because the 
coefficient C n in eq.(45) behaves as n\ at large n. 

The first problem is mainly technical and has been solved by the BFKL equation. 
However this equation carries with two problems, namely, the next order corrections have 
yet not been found and the solution of the BFKL equation is influenced by confinement 
region. Moreover, such confinement corrections have not been localized in something like 
the matrix element of some operator. 

The second problem is much more complicated one and it can be solved only if we 
understand better the nonperturbative origin (renormalons, instantons,etc) of behaviour 
of coefficients C n . We are only in the beginning of the systematic study of the nonper- 
turbative effects in high energy problem ( see ref. jffSH ). 

The origin of the shadowing corrections is also nonperturbative one, but it is new kind 
of nonperturbative phenomena which closely related to the high density of the partonic 
system. At the moment we can control them theoretically only in the transition region 
between pQCD and hdQCD in Fig. 19. 

Coming back to the list of our hopes in section 9.5 we can conclude that we have 
understood better the kinematic region in which we can trust pQCD. 
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16 Nearest Future = Theory in the transition region. 



I firmly believe that in the nearest future : 
1. 

a theory for the transition region between pQCD and hdQCD ( see Fig. 19) based on 
the BFKL Pomeron and Pomeron interactions will be completed. 

2. 

we will understand better the manifestation of the BFKL Pomeron and SC (shadowing 
corrections) and coming experiments from HERA and Tevatron will help us to understand 
better in which kinematic region we can restrict ourselves by new Reggeon Calculus based 
on the BFKL Pomeron. 

3. 

we will find the values of all nonperturbative parameters that we need (like R in the 
GLR equation) with help of the nonperturbative QCD methods such as lattice calculation 
or QCD sum rules. 

The Mueller approach as well as the deep understanding of the correct degrees of 
freedom allows us to start to built the effective theory for high energy interaction. It is 
interesting to emphasize that this effective theory will be a theory of interaction of the 
dimensional objects ( colour dipoles). Therefore, at high energy we have such a system of 
partons (particle) in which we can study the theory of interaction of the dimension degrees 
of freedom. It seems very attractive, because we can check our theoretical imagination 
with experimental data, while the string gang who is doing basically the same trying to 
understand how can we build the relativistic theory for dimension degrees of freedom still 
have to rely only on theory without any experimental support. 

17 Toward the effective theory. 

I hope that I have convinced you that the main direction of our movement is toward the 
effective theory based on nonperturbative QCD approach, because the principle problems 
of high density partonic system cannot be solved without understanding of npQCD. We 
hope that in the future we will find: 

1. a selfconsistent effective theory for a high density system in QCD; 

2. the analytic solution of hdQCD, using the smallness of the coupling QCD constant 
in this kinematic region; 

3. a new collective phenomena in hdQCD. 

At the present we can see two harbingers of this shining future: 
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17.1 Effective Lagrangian for high energy QCD. 



The first attempt to develop nonperturbative approach to low xb problems was to write 
down the effective Lagrangian for the region of small xb- Intrinsically we assume that 
such a Lagrangian should be simpler that the QCD one and allows one to apply some 
direct numerical procedure (lattice calculation for example) to calculate the scattering 
amplitude with this Lagrangian. It should be stressed the attempts to calculate the 
amplitude with full QCD Lagrangian have failed by now. At the moment we have two 
effective theories on the market: one was proposed by Lipatov ^3] which looks not much 
simpler that full QCD but it certainly incorporates all results of perturbative calculations, 



and the second was suggested by Verlinde and Verlinde p0[, which is much simpler and 



is suited for lattice -like calculation but it has not been checked how well this effective 
theory describes the perturbative results. Moreover there is some indication |)TJ that this 
theory cannot describe the virtual correction in the BFKL Pomeron. 



17.2 Thermodynamics of high density QCD. 

I firmly believe that we need to write down the correct kinetic equation for high density 
QCD. Such an approach has certainly at least one big advantage: the smooth matching 
with the GLR equation. Unfortunately we have not yet understood how to write such 
an equation in our nonequilibrium situation. However we have understood better the 
physical meaning of the new typical momentum in our parton cascade ((|pt|) in section 
11.2 ). It turns out that this momentum is the Landau - Pomeranchuk momentum for 
our parton medium p4|. Thus the gluon emission with transverse momentum less that 



(\pt\) is small due to destructive interference between emission before and after collision 
of the parton with other partons in the medium. The experience with the calculation of 
the anomalous dimension of high twist gluonic operators also give us understanding why 
for bosonic degrees of freedom such as a gluon there could be saturation of gluon density. 
Indeed, our cascade is rather a one dimensional one and for such a system the direction 
of motion plays the role of spin for fermions. 



18 Conclusions. 

We can perhaps summarize our discussion with the following statements. 

1. The new kinematic region of high density QCD is wide open for theoretical and 
experimental investigation. 
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2. Perturbative QCD methods have been developed to explore the transition region 
between the region routinely analyzed by perturbative calculations for hard processes and 
the hdQCD region. 

3. Between the kinematic region of soft processes (key words: unitarity, analyticity, 
Reggeons, Reggeon Calculus, hadronization models, etc) and the kinematical region of 
hard processes (key words: quarks, gluons, parton distributions, asymptotic freedom, 
pQCD, evolution equations, etc) is the region of semi-hard processes, where the cross 
sections can be as large as the geometrical size of the hadron but with a typical transverse 
momentum large enough that methods of pQCD can be applied. 

4. New probes have been developed for this high-density QCD region. Although we 
still lack some theoretical understanding and cannot provide reliable estimates for these 
probes, we firmly believe that they provide a deeper insight into QCD and could lead to 
new experimental discoveries. 

Phase Transition. 

We can represent the predictions of the previous section as a phase transition between 
different forms of hadronic (partonic) matter. Our prediction can then be summarized as 
the possibility of a phase transition from a quark - gluon ideal gas to a gluon liquid (see 
Fig. 36). 

Indeed, in the pQCD region (see Fig. 19) we can regard a partonic system as an ideal gas 
of quarks and gluons because we can neglect their interactions (only emission is essential 
here) . 

When the density of partons increases we approach the region where we must consider 
the interactions between partons in a simple approximation and in a situation where this 
interaction is not too big. This then is a Van der Waals gas of quarks and gluons. Here 
we have a theory based on a nonlinear evolution equation that is the analogue of the Van 
der Waals equation for a parton system. 

For even larger values of the density, the smallness of the coupling constant results 
a small correlation between partons. We call this system of partons a gluon liquid. We 
recall that the density is high primarily because of emission of gluons and quarks inside 
this system, which still can be treated in a perturbative way. 

The border between the gluon liquid and the Van der Waals gas is the critical line we 
discussed in Section 12.3. The trajectory of the GLAP evolution equation that touches 
the critical line separates the Van der Waals gas from the ideal gas. The asymptotic of 
these borders are also calculable. 

In conclusion we want also to recommend you the list of the reviews on the subject 
[ p^j which, we hope, you will be able to understand after these lectures. 
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Figure Caption: 

Fig.l: The exchange of the resonance with spin j. 

Fig. 2: The exchange of the Reggeon. 

Fig. 3: The Reggeon trajectories. 

Fig. 4: Resonance - Reggeon Duality. 

Fig. 5: The Pomeron structure in the Veneziano model. 

Fig. 6: The structure of the parton cascade. 

Fig. 7: The random walk of parton in transverse plane. 

Fig. 8: The single diffraction dissociation process in the Pomeron approach 

Fig. 9: The triple Pomeron diagram. 

Fig. 10: An example of enchanced Pomeron diagrams. 

Fig. 11 The AGK cutting rules for two Pomeron exchange. 

Fig. 12: Hadron -dueteron interaction. 

Fig. 13: The process with multiplicity 2 < n N > (a) and < > (b). 

Fig. 14: The v Pomeron exchange diagram with jx cut Pomeron. 

Fig. 15: Pomeron diagrams for the total,elastic and diffraction dissociation cross sections. 

Fig. 16: Pomeron diagrams for inclusive cross section (a) and rapidity correlations (b). 

Fig. 17: One and two parton shower processes. 
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Fig. 18: Single diffraction dissociation in Reggeon approach. 

Fig. 19: The map of QCD.p is the parton (gluon ) density 
and r is the distances resolved in an experiment. 

Fig. 20: The "Pomeron" Calculus in QCD. 

Fig.21: The Born Approximation (BA) of pQCD. 

Fig. 22: The next to BA of pQCD for quark-quark scattering: 
emission of one extra gluon. 

Fig. 23: The next to BA of pQCD for quark-quark scattering: 
a§ correction to elastic amplitude. 

Fig. 24: The next to BA of pQCD for quark -quark scattering: 
gauge invariance trick. 

Fig. 25: The next to BA of pQCD for quark -quark scattering: 
the resulting answer for one extra gluon emission. 

Fig. 26: The relation between colour coefficients. 

Fig. 27: The BFKL equation. 

Fig. 28: The basic branching process in DIS. 

Fig. 29: The parton distribution in the transverse plane. 

Fig. 30: Structure of the parton cascades at low x and 
the coherence in "ladder" diagrams. 
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Fig. 31: The random walk picture in Ink 2 . 

Fig. 32: The behaviour of xG(x, Q 2 ) versus Q 2 at fixed x. 

Fig. 33: Trajectory of the nonlinear GLR equation. 

Fig. 34: The Mueller - Navalet process. 

Fig. 35: Large Rapidity Gap process. 

Fig. 36: The phase diagram for hadron - hadron collisions for zero rapidity 
in cm frame. W is the absorption probability defined in eq. (103). 
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